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A MODEL OF THE EFFECTS OF IRRADIATION ON BACTERIAL 
ABILITY TO SUPPORT PHAGE GROWTH* 


BE. T. EPSTEIN 
BRANDEIS UNIVERSITY 
WALTHAM, MASSACHUSETTS 


It is pointed out that within a phage-infected bacterium there may be 
several types of entities which can permit such a cell to produce plaques 
when plated before infection. From the equation giving the survival of 
the cell as a plaque-producing entity it is shown that actual survival 
curves will depend markedly on the relative radiosensitivity of the 
various entities. Survival curves in the literature are examined and 
interpreted as revealing the existence of several phage-reproducing 
factors, in unactivated form, in the cell before infection. The following 
model is shown to be consistent with all the data examined. After in- 
fection, the phage and bacterium undergo modifications lasting about one- 
third of the latent period. After this time replication begins with the 
aid of the phage-reproducing factors which are successively activated 
until all are in production; this takes until the end of the first half of the 
latent period. After this time the pattern of survival curves depends on 
the relative radiosensitivity of the phage-reproducing factors and the 
phages which begin to appear at about two-thirds of the latent period. 


S, E. Luria and R. Latarjet (1947) initiated the radiobiological 
study of intracellular phage synthesis by irradiating infected 
bacteria at various times in the latent period and plating the bac- 
teria before the end of the latent period to determine the survival 
of the bacteria as phage-producing entities. In their analysis of 
the results, these authors pointed out the many deviations from 
the perhaps-to-be-expected multitarget survival curves based on the 
notion that starting from one infecting phage there should be some 
time at which there are two, three, etc., progeny as the infecting 
s its replication. R. Latarjet (1948), S. Benzer (1952), 


hage effect 
ie Jacob (1953) did similar studies and found 
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at least three different patterns of radiosensitivity changes, lead- 
ing to the notion that there may be three different modes of phage 
reproduction. Because of these unexpected types of results, 
radiobiology has not yet been able to contribute to the elucidation 
of the sequence of events in virus reproduction. The purpose of 
this paper is to show that a slightly more general approach to the 
interpretations of the results of these investigations reveals a 
hitherto unsuspected pattern in the data, suggesting a new mech- 
anism in the reproduction. 

First we write the expression for the survival of an infected 
cell as a phage-producing entity. Within such a cell there may 
be simultaneously mature phages, embryonic phages which have 
only to be assembled, phage synthesizing mechanisms, and per- 
haps other unsuspected entities which can, by themselves, lead to 
phage activity. If p, is the probability that the ith type of entity 
survives a given dose of radiation D, (1 - p,)"t is the probability 
that all n; copies are inactivated by the irradiation. If there is 
more than one type of entity which can lead to phage activity, 
the probability that all such types are inactivated is the product 
of several such terms, and the probability that at least one of 
these entities survives is one minus that product, Also, p; may 
be written as exp(—D/D ;) where D,; is the 37% survival dose for the 
ith type. Thus, finally, the surviving fraction (N/N,) of phage- 
producing cells is 


N/N, =o (te =D aP\es (l-e eGR EAR. ieee ~D/Dkyrk, 


The important features of this relation are that for high doses 
the asymptotic slope of the curve log N/Ny vs D is that of the 
most radioresistant entity and extrapolation of this slope back to 
zero dose gives as intercept the number n of these entities, These 
features follow directly from the equation by letting D become very 
large. Having determined the asymptotic slope and the asymptotic 
multiplicity from a survival curve, one could, in principle, deter- 
mine the number and radiosensitivity of the next most resistant 
entity by subtracting from the survival curve the survival due to 
the most resistant type of entity. In practice, however, the data 
are not ordinarily accurate enough to permit more than a qualitative 
conclusion about the presence of other types of entities. 

This analysis will be applied first to the survival of the T2-F, 
coli B system studied by Luria and Latarjet (Joc. cit.) and by 
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Benzer (oc. cit.) using ultraviolet light (UV) and by Latarjet (Joc. 
cit.) using xerays. The x-ray data are given in Figure 1: in this 
figure our additions are the dashed extrapolatory lines and the 
circles. From the figure it is seen that the multiplicity (n,) is six 
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FIGURE 1. The fractional survival of T2-infected E.coli B at varying 
doses of x-rays. Data are from Latarjet (1948). The number on each 
curve is the time in minutes since infection began. The circles are points 


of a six-target curve calculated to have the same asymptotic slope as the 


curve near which the circles lie. 


ntire second half of the latent period. (From other 


throughout the e 
data Filion by Latarjet it is found that n, increases from unity 


starting at about % the 21 minute latent period and reaches two at 
9 minutes). It is well-known (Doermann, 1952) that there are ai 
mature intracellular phages at the start of the second half of the 
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latent period and that there are 30 or 40 at 18 minutes. Thus it 
follows that the entities revealed by the irradiation cannot be 


mature phages. ; 
Because of the technical difficulties involved in such experi- 


ments, the number six should be construed as meaning, say, a 
number between 3 and 10. However, the number six is also found 
as the extrapolate multiplicity in a similar analysis of the UV 
studies of the same system by Luria and Latarjet (loc. cit.) and 
by Benzer (loc. citt.). 


PLAQUE-FORMING ADSORPTIONS 


10 20 30 40 
MINS UV 


FIGURE 2. The fraction of adsorbed T2r phages which produce plaques 
on E.coli B irradiated with ultraviolet light for the times shown. The 
solid line is a six-target curve chosen for best fit to the data. N/N» = 
L— (1, =e" P/P1) 4s ee be Dek minke LXE 


The next question considered is the nature of the six entities, 
In this connection we ask if the entities exist within cells before 
infection. This can be answered by determining whether a six- 
target curve is obtained for the “‘capacity’’® of the cells. The 
Capacity is the ability of irradiated cells to yield phage progeny 
if infected after irradiation. A theoretical six-target curve for the 
T2-E.coli B system is given in Figure 2 as the solid line, as- 
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suming only one type of entity exists. A reasonable agreement 
with the data is evident. Thus it is likely that the multiplicity is 
characteristic of some cellular structures or factors; they will be 
called phage-reproducing factors (prf’s). A possible interpretation 
of the prf’s is that they are the sites of phage replication. 

The model being developed thus predicts that the survival curves 
should increase from unit multiplicity after the first prf has been 
activated by the infecting phage. After they are all activated, the 
survival curves’ asymptotic slopes and multiplicities should 
depend on the relative radiosensitivity of prf and maturing phage. 
If the phage is more radiosensitive, the survival curve should 
remain essentially that of the prf’s. If the prf is more radiosensi- 
tive, the asymptotic slope and multiplicity should be that of the 
phages, and this change should become evident when the first 
phage appears, starting at about % of the latent period. 

There is one other system for which both survival curves and 
the capacity curve are in the literature: the p8-P.pyocyanea sys- 
tem studied by Benzer and Jacob (loc. cit.) In this case the 
phages are less radiosensitive than the sites, and the sequence of 
slopes and multiplicities is just as predicted in the preceding 
paragraph for such a case. The authors determined the intracel- 
lular phages by artificial lysis of the bacteria and pointed out the 
correlation of their appearance with the changing character of the 
survival curves. The capacity multiplicity is two. 

Studies on two other systems have been reported in the literature 
which has been examined. Benzer’s study (loc. cit.) of the T7-E. 
coli B system with UV shows the features characteristic of the 
case of the phage being more radiosensitive than the prf. The 
extrapolate multiplicity is six. Benzer and Jacob (loc. cit.) stud- 
ied the phage lambda-E. coli K12 system. It is similar to the T2 
case, but the results and analysis are not as cleancut for the data 
stop at about the time when the first mature phages should appear. 
In neither of these two cases was the capacity determined; so 4 
critical test of that aspect of the model cannot yet be made. 


Discussion 
It appears reasonable to conclude that it is not valid to attempt 


to interpret the radiation survival curves of virus-infected cells 
solely on the basis of the existence of one type of phage-preducing 


entity, for the curves either show directly the existence of more 
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than one type of entity or else reveal an entity which clearly 
cannot be a phage, 

The model proposed is so general as to be only slightly re- 
strictive, but it shows a hitherto unknown pattern of development 
in the survival curves. The observation of this pattern makes it 
unnecessary to suppose three markedly different modes of phage 
synthesis as might otherwise be proposed on the basis of the 
markedly different sequences of survival curves. 

Finally it is noted that the prf model affords a mechanism for 
post-infection interference phenomena for, if all the prf’s are acti-e 
vated by one phage strain before a second strain can reach the 
prf’s, the second strain should effectively be excluded even though 
it might otherwise develop within the particular cell. 
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CONTRIBUTIONS TO THE TEEORY OF ACTIVE TRANSPORT 


CLIFFORD S. PATLAK 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO* 


A system consisting of two solutions separated by a membrane may be 
in one of four possible states: (1) transient, (2) steady, (3) equilibrium, 
or (4) pseudo-equilibrium. The latter state denotes that in the solutions 
the net flow of all components is zero but at least one of the components 
is not in thermodynamic equilibrium. Transient and steady-state systems 
may or may not have active transport. Thus only systems in either equi- 
librium or pseudo-equilibrium are considered in this paper, since the 
former indicates that there is no active transport, whereas in the latter 
case there always is active transport. This simplifies the problem of 
finding whether a system does or does not have an active transport 
mechanism, since it is frequently fairly easy to determine experimentally 
whether a system is in equilibrium or pseudo-equilibrium. 

The assumption that electric neutrality exists within very thin mem- 
branes is shown not to be valid. However, electric neutrality does exist 
in the solutions in a system in a pseudo-equilibrium state with fixed 
charges and impermeative ions. It is then shown how the presence and 
sign of an electric potential may be found by use of electroneutrality. 

The mechanism of active transport may be due to a general force acting 
on all particles of a particular component or to an individual force acting 
on the individual particles of a particular component. A general solvent 
flow or a diffusion drag force illustrates the first mechanism while the 
second is accomplished by either a carrier or a Maxwell Demon. The 
general type of active transport has been extensively treated in the 
literature, while the individual type has not been treated in a generalized 
form. Therefore, the individual type of active transport is discussed at 
length, and a simple illustrative model is intensively analyzed. Follow- 
ing this, there is a discussion of the Maxwell Demon and some models of 


it are presented. 


Introduction: 

It has been observed in many biological situations that the net 
transport of materiai from one region to another cannot be explained 
as due to a passive movement following thermodynamic laws. The 


*This research was supported by the United States Air Force through 
the Air Force Office of Scientific Research of the Air Research and 
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term ‘‘active transport’? has been applied to such cases and will be 
discussed in this paper after the introduction of some necessary 
definitions and nomenclature. 

We shall be concerned solely with systems of uniform temperature 
and-pressure which are composed of two different solutions sepa- 
rated by a membrane and in which the properties of the two so- 
lutions are known. We shall use the term molecule to signify either 
a charged or uncharged particle. The solutions then contain various 
molecules which are freely diffusible between the two solutions, 
molecules whose flows are linked, molecules which are fixed (such 
as arise from binding on proteins), and molecules which are not 
fixed but cannot diffuse between the solutions. By diffusing be- 
tween the two solutions we mean the ability to passively traverse 
the membrane (e.g., by actual diffusion or via some type of non- 
energetic carrier), By. linked flow we mean that the movement, 
across the membrane in either direction, of one or a group of mole- 
cules is linked together. Both of these terms will be extensively 
examined in the discussion of permeability at the end of this section. 
The membrane will be permeable to the diffusible molecules and 
may possess properties which allow it to actively transport some of 
the molecules (both diffusible and not). We define an empermeative 
molecule to be one which cannot cross the membrane and also is not 
actively transported across it. We define a permeative molecule to 
be one which is not impermeative (i.e., can traverse the membrane 
by some mechanism) and does not have a linked flow. A component 
of the system shall be a permeative molecule. 

An individual component of a system may be in one of four states, 
determined by j (X,2), which is the net flux of the component under 
consideration across a unit area at time ¢ and position X, 

1, The transient state. In this case j(X,) is a non-constant 
function of time. This situation arises, for example, in a nerve or 
muscle fiber when it is excited, in a receptor organ (such as the 
retina) when stimuli are present, etc. 

2. The steady state. In this case, j(X,z) is a function of X alone 
(i.e., for each X it is a constant with respect to time) not equal to 
zero at all points of the system. This situation arises, for example, 
in almost all diffusion flows across a membrane after a sufficient 
time has passed for the transient state to disappear. (Here j(X,d) 
is not equal to zero in the region of the membrane but may equal 
zero in the two ‘‘reservoir” regions on either side of the membrane.) 
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3. The equilibrium state. In this case j(X,¢) is equal to zero for 
all X, and the ratio of the concentration of the component at X, to 
that at X, is determined completely by the characteristics of the 
points X, and X,. That is, the concentration ratio is determined 
by the equality of the partial molar free energy of the components 
at the two points. As is well known, the partial molar free energy 
at a point is determined solely by the region immediately surround- 
ing the point. This situation arises, for example, in a Donnan 
equilibrium. 

4, The pseudo-equilibrium state. In this case j(X,¢) is zero, as 
in the equilibrium state. However, the concentration ratio of the 
component between two points is not determined solely by the 
region immediately surrounding the two points but is also a function 
of the characteristics of the remainder of the system. This situation 
arises in the case of an active transport where the component is 
being actively transported across a membrane at the same rate as 
it is diffusing back across the membrane so that the net flow across 
any surface is zero. This situation is probably found in many 
organs in vivo that have a resting potential, such as peripheral 
nerve fiber, rods and cones in the retina, etc. It appears that the 
term pseudo- equilibrium used here is equivalent to the term ‘‘dynamic 
equilibrium’’ used by R. Hdber (1945, p. 244). 


The above definitions concern the state of an individual com- 
ponent. We define the state of the system to be the lowest order 
state found among the individual components of the system. To do 
this, we use the conventional approach of ordering the states of 
the individual components in the following ascending order: (1) tran- 
sient, (2) steady, (3) pseudo-equilibrium, (4) equilibrium. Thus, 
the resting nerve fiber én vzvo may be in a pseudo-equilibrium state; 
the same fiber in vitro is in the steady-state (H. Grundfest 1955), 
while the active nerve fiber is in the transient state. We shall 
later in this paper talk of putting the system into at least a pseudo- 
equilibrium state. By this we mean that we put the system into a 
state such that j(X,¢) will equal zero everywhere for all components, 
and therefore it may turn out to be either in pseudo-equilibrium or 
in true equilibrium. 

A major problem of biology has been the attempt to decide whether 
a component in a given system is undergoing active transport. 
When the system is in the steady state, exceedingly difficult prob- 
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lems arise that have not been solved in a satisfactory manner. 
T, Rosenberg (1948, 1954) has stated that the fact that the net flow 
of a material is from a region of low electrochemical potential to 
that of a higher one is evidence of active transport. However, as 
Rosenberg acknowledges, this is not general enough, as will be 
Shown in the following four examples. 

As pointed out by D. Spanner (1954), E. Conway (1954), and 
H. Ussing (1954), the component may be actively transported in the 
direction of high to low electrochemical potential such that its net 
flow is greater than if it were passively diffusing. We may extend 
this and note that the component may also be actively transported 
from low to high electrochemical potential at a rate lower than 
that at which it is diffusing in the opposite direction. Then its net 
flow is less than if it were only passively diffusing. In both cases 
the net flow will be from a region of high to a region of low electro- 
chemical potential. 

Another example in which Rosenberg’s criterion would not work 
is the case of a closed system with no production or consumption 
and with two components, the first being present in much higher 
concentration than the second, Let the distribution of the two 
components be such that where the electrochemical potential of 
the first component is high, the electrochemical potential of the 
second is low, and vice versa. Let there be a large diffusion 
drag force between the two components, and let the membrane per- 
meability to the two components be approximately equal. Then the 
second component would appear to be diffusing from a region of low 
to that of a high electrochemical potential. If the regions are large 
enough and the rates of flow small enough so that for a sufficient 
length of time the system appears to be in a steady state, then by 
Rosenberg’s criterion there should be active transport. We will 
assume that there must be an expenditure of energy in order for 
active transport to occur. Although the criterion for active trans- 
port can be made quite explicit assuming energy need not be ex- 
pended for the active transport (D. Spanner, Joc. cit.), the term 
active transport would lose much value and usefulness, In the 
previous example, since energy is not expended and the movement 
1S passive, we see that there is no active transport. 

A third example is a modification of the above model, in which 
the components travel across the membrane by means of a non- 
energetic carrier. To specify this model we assume that (a) we 
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have a closed system with no production or consumption; (v) the 
characteristics of the two components with respect to the carrier 
are identical; (c) the carrier is present in such a low concentration 
that it is always saturated but it can traverse the membrane without 
carrying any component; (d) the components are neutral so that the 
electrochemical potential is measured solely by the concentrations; 
(e) the carrier ‘‘loads’® and ‘‘unloads’’ only at the boundaries; and 
(f) on one side of the membrane the second component is present in 
much lower concentration than the first component, whereas on the 
other side of the membrane it is present at a higher concentration 
than the first component, although at a lower concentration than its 
value on the first side. In this case, the second component will 
not be carried across from its high concentration to its low concen- 
tration due to the ‘‘swamping’’ effect of the first component. How- 
ever, it will be able to go across in the opposite direction, since it 
is present there at a higher concentration than the first component. 
Thus, there will be a net transfer from low to high electrochemical 
potential and, by the same arguments used in the previous example, 
we have a violation of Rosenberg’s criterion. 

A fourth example is anomalous osmosis (R. Héber, loc. cit., 
pp. 624 ff.), the explanation of which is fundamentally similar to 
that of the second example. In the last three examples the apparent 
active transport fundamentally arose from the fact that the system 
is not in a ¢rue steady-state but only in a semi-steady-state. That 
is, if a sufficiently long time is allowed to pass, it would be ob- 
served in the latter three examples that the direction of flow would 
have changed and would go in the direction expected from thermo- 
dynamic considerations. However, this factor cannot be dismissed 
lightly but must be considered in some manner, since experiments 
can only run for a finite time due to the limited time that the @X- 
perimental material will last and the limited time that the experi- 
menter can devote to the experiment. 

We do not have a better criterion than Rosenberg’s to offer for 
the: case of the steady-state. However, we do have a criterion to 
apply to systems in true or pseudo-equilibrium. lf the system is in 
at least a pseudo-equilibrium state, a determination of whether any 
component is in pseudo-equilibrium or whether all compenents are 
in true equilibrium can be used to determine the existence or non- 
existence of active transport. If any component 1S 1n pseudo- 
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equilibrium, active transport must exist; whereas if all components 
are in true equilibrium, active transport is not present. 

Some biological material (e.g., frog-skin; cf. E. Huf, 1955) can 
be placed in at least a pseudo-equilibrium state. However, there 
are other materials (e.g., nerve fiber) with which it is more difficult 
to work, Although it may be conceptually possible to place them in 
at least a pseudo-equilibrium state, they cannot be so modified due 
to the tremendous experimental difficulty attached to the procedure. 
Also, to do so would in many cases necessitate changing the con- 
centration of the components to such a large extent that the system 
would no longer be in a normal physiological state. A secreting 
tissue (e.g., mucosa of the stomach) is an extreme example of a 
tissue that would be drastically disturbed by any attempt to place 
it in at least a pseudo-equilibrium state. 

In many cases it is comparatively easy to determine, with the use 
of tracers and other means, whether the system is in at least a 
pseudo-equilibrium state (as contrasted to a steady or transient 
State). However, knowing the system is in at least pseudo-equi- 
librium, it is more difficult to determine whether any component of 
the system is in true or in pseudo-equilibrium. If the solution on 
both sides of the membrane can be controlled, as in the frog skin 
experiments, and if it is known whether or not a molecule is perme- 
ative, this problem does not arise. However, if one of the two 
Solutions cannot be controlled (such as the intracellular fluid in 
muscle), it then becomes more of a problem to decide whether the 
system is in true or pseudo-equilibrium. [This problem arose in 
the debate concerning the **fixed-charge’’ theory of G. Ling (1952). 
It also arose in the attempt to decide whether there was a true or a 
pseudo-equilibrium between the aqueous humour and the plasma of 
the eye (H. Davson, W. Duke-Elder, and D. Maurice, 1949).] The 
impossibility of differentiating between equilibrium and pseudo- 
equilibrium in certain cases is the basis of a theory of apparent 
active transport proposed by A. Bierman (1953), This general 
problem will be discussed further in the section on types of active 
transport. For simplicity, throughout the remainder of this paper 
we shall assume that all properties of the solutions are known. 

It should be emphasized that, even if no active transport is found 
as a result of our analysis (i.e., the system is in the true equilibrium 
state), this of course does not prove that the system does not pos- 
Sess an active transport mechanism when it is in a different (and 
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lower) state, such as the steady state. For example, consider a 
System which normally has some material flowing through one of 
the regions. Let this material combine with one of the components 
of the system. Then let the combined form be carried into or across 
the membrane. At this point let the combined form be broken up 
and the material consumed. Thus, we have an active transport 
mechanism which would no¢ appear in a true or pseudo-equilibrium 
state because the experimental set-up would not allow for the 
introduction of this ‘‘carrier’? material. Further examples are ob- 
tained as modifications of the models presented in the previous 
discussion of the Rosenberg criterion. Assume that material ap- 
pears in a solution by production within that solution or flow from 
another solution. Assume further that the same material disappears 
in the other solution by consumption within that solution or flow to 
another solution. In this manner the necessary conditions for the 
two examples to exist in a true steady-state are met. Then we 
have a type of active transport which would not appear in an ex- 
periment in which there is at least a pseudo-equilibrium state. The 
mechanism of the second example in the discussion of the Rosenberg 
criterion has been extensively discussed by J. Hearon (1950, 1953). 
In all cases, there would be an energy expenditure involved in the 
inflow-outflow or production-consumption mechanisms. Still another 
example would be that in which the temperature of the two regions 
differs (D. Spanner, loc. cit.). 

We can limit the term ‘‘active transport’? to apply only to those 
cases in which the mechanism is contained completely or partially 
within the membrane. Then the above examples would not have 
been examples of active transport. However, we have chosen not 
to do this. In any case it would make no difference for the pseudo- 
equilibrium state. As these examples illustrate, if there is any 
active transport in the pseudo-equilibrium state, the mechanism 
would have to reside within the membrane. 

Whether or not active transport is present in the steady-state 
system can in certain cases be determined best by the pseudo- 
oach. ‘That is, if there is an energy expenditure 
t in the steady state, but if the system 
tends to go with increasing time to a true equilibrium system, it is 
perhaps best not to call this active transport. An anatenes. of this 
can be found in the earlier example of the non-energetic carrier, in 
which it is entirely possible that energy may have been dissipated 


equilibrium appr 
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by the system. This is related to the discussion of ‘facilitated 
diffusion’? by J. Danielli (1954). However, since this mechanism 
would actually yield a true equilibrium, we feel that this should 
not be called active transport (as opposed to Danielli, loc. cit.) 
even though energy is being dissipated. 

We will now discuss criteria useful in determining whether or 
not a molecule is permeative. We might say a molecule is permeative 
if it can traverse the membrane (as evidenced by a tracer of the 
molecule appearing in the solution on one side of the membrane 
after it has been introduced into the solution on the other side) and 
impermeative if it cannot traverse the membrane. However, this 
simple approach leads to many difficulties. After discussing three 
general examples to illustrate these difficulties, we shall restate 
our definition of permeativity. This discussion will show the ne- 
cessity for our definition. 

All three of the examples have certain common features. One of 
these is that the molecule can traverse the membrane, as evidenced 
by tracer studies, but not by simple passive diffusion. Also, when 
the system as a whole goes to at least a pseudo-equilibrium state, 
the molecules crossing the membrane in the three examples do not 
go to equilibrium. We have said that when a system is in at least 
a pseudo-equilibrium state, if all of the permeative ions are in 
equilibrium the system has no active transport whereas if they 
are not all in equilibrium the system does have active transport. 
However, since there is no active transport in these examples and 
since the molecules are no¢ in equilibrium, the above criterion 
would be violated. For this reason we would wish to say that the 
molecules in these examples are not permeative. Hence, we must 
find some factor which is common to all of these examples (besides 
the fact that they do not go to equilibrium). This factor, that the 
flows are linked, will be discussed after the three examples, 

The first example is the case of simple exchange diffusion as 
defined by H. Ussing (1949c). For this case a molecule, Q, at- 
taches itself to a carrier on one side of the membrane (the internal 
Solution side), is carried across the membrane, and is then unloaded 
into the solution on the other side of the membrane (the external 
Solution side), Then, in order for the carrier molecule to return to 
the internal solution side of the membrane, it must combine with 
another Q and carry it back across. Similarly, for the carrier on 
the internal solution side of the membrane to go to the external 
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solution side of the membrane, it must also combine with a mole- 
cule of Q. If we assume that Q cannot cross the membrane by any 
other method than the above, then, except for the transient state, 
for each molecule of Q that is carried from one side of the mem- 
brane to the other, another molecule of Q must be carried in the Oop- 
posite direction. We thus see that there cannot be any net trans- 
port of the molecule across the membrane regardless of the mole- 
cular distribution in the solutions on either side of the membrane. 
However, if a tracer is introduced into one of the solutions, it 
will traverse the membrane and appear on the other side. For any 
initial distribution, the final distribution of Q will be the same as 
the initial distribution. Hence, when the system is in at least a 
pseudo-equilibrium state, the distribution of Q will not be deter- 
mined by equating the electrochemical potential of Q on the two 
sides of the membrane. 

The second example which we shall consider is a generalization 
of the above example. In this case, instead of the carrier trans- 
porting only one type of molecule, it may transport many types of 
molecules, For definiteness, we will assume that it can transport 
three types of molecules, Q, R&, and S, and that the carrier may 
traverse the membrane only when carrying one of them. Further, we 
will assume that Q, R, and S cannot cross the membrane by any 
other method. Thus, considering Q, R, and S as a unit, there is no 
net transport of the three types across the membrane. However, if 
tracers of Q, R, and S are introduced into one of the solutions, 
each will traverse the membrane and appear on the other side. 
When the system is in at least a pseudo-equilibrium state, the 
distribution of the Q, R, or S alone on either side of the membrane 
will of course be different than their initial distribution. However, 
the sum of the Q, R, and S on either side of the membrane will be 
the same as the sum for the initial distribution. Therefore, the 
distribution for the Q, R, and S will not be determined by equating 
the electrochemical potentials on either side of the membrane. If 
we alter our assumptions by permitting one of the molecules to 
traverse the membrane by some other means, then there can now be 
a net transport of Q, R, and S acting as a) unit. In this case, the 
three types of molecules can go to an equilibrium state for each of 
the molecules individually. That is, even though @ and # can only 


traverse the membrane via this carrier, the fact that S can also 
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traverse the membrane via a different method now allows all three 
molecules to go to an equilibrium state. 

The third example we shall consider was proposed by M. Jacobs 
(1940) and W. Osterhout (1940) and treated under the term “‘partial 
equilibrium’? by T. Rosenberg (1948). We shall consider a specific 
case under this example involving three molecules, Q, R, and S, 
This mechanism assumes that Q, FR, and S cannot cross the membrane 
separately, but their chemical combination, QRS, can. Thus, for 
each molecule of Q that crosses the membrane, in either direction, 
a molecule of # and a molecule of §S must also cross in the same 
direction. Therefore, tracers of Q, R, and S will traverse the 
membrane. When the system is in at least a pseudo-equilibrium 
state, QRS will be in equilibrium. However, the concentration of 
Q, R, and S separately will not be found by equating the electro- 
chemical potentials of each of the molecules in the two solutions. 
If we modify our assumptions so that S is now able to traverse the 
membrane by some other method, then in at least a pseudo-equilibrium 
state S would be in equilibrium. However, Q and RF still would 
not be in equilibrium. This is in contrast to the result of the 
second example on a generalized exchange-diffusion, where S was 
assumed able to traverse the membrane by some other method. 

A common factor for all of the above examples is that the flows 
of the molecules are no¢ independent but are linked. We shall use 
this particular factor because, among other reasons, it may, in 
principle, be experimentally tested in a straightforward manner. 

In order to determine whether or not the flows of several mole- 
cules are independent we will use the following procedure, after 
introducing some nomenclature. Let the subscript 7 denote the /th 
molecule, let j and ; denote the flow per unit time and area across 
the membrane from the internal to the external solution and the 
flow from the external to the internal solution respectively, and let 
v be a constant. 

For at least the steady state the values of i; and a for all of the 
molecules may be found by means of tracers. Consider the equation 


>> ee 
Li + 25%), = 9, (1) 
I 1 


: Ti > ie 
for all possible values of the j, and j in which both summations 
are over all the molecules. Then, altering the concentrations but 
still having the system remain in a physiological state, we find 
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> + 
what the new values are of j, and j,. Doing this a sufficient number 
of times, we see whether we can find values of v, and y not all 
equal to zero and not in the obvious combination that can always 
be found to express electric neutrality, such that equation (1) is 
satisfied for al? the different possible concentration distributions. 
If constants v, and y can be found to satisfy equation (1), then all 
of the molecules whose coefficients are not equal to zero are linked. 
(Mathematically, we may say that if equation (1) is satisfied for 
any set of /’s, then the j,’s of that set are not linearly independent). 
Let us evaluate, by the criterion of equation (1), the three previous 
examples to see if they have linked flows. 

For the first example of simple exchange diffusion, equation (1) 

would be 

= 

} ag fo = 0. 
For the second example of a more general exchange diffusion, in 
which Q, #, and S can only cross the membrane via the carrier, 
equation (1) would be 


Pee ee ee 
Gio + ip t is)— Ug + Jp +45) = 9. 


Note that if S is also allowed to traverse the membrane via a different 
route, the above equation is no longer valid. Therefore, the mole- 
cules would no longer be linked, which is the same result arrived 
at above. 

For the third example of chemical combination, in which Q, RF, 
and S can only cross the membrane as the compound QS, equation 
(1) would be 


If S is also allowed to traverse the membrane via a different route, 
this will not affect the equation involving @ and #&. Therefore, 
they will remain linked, as we had concluded above. 

Obviously, the equations for all three of the above examples 
would be valid for add concentration distributions of Q, &, and S in 
both solutions. Thus we can say that the flows of Q, R, and S are 
linked, a result we had arrived at before by a different argument. 
However, we see that the criterion of equation (1) is the result of a 


282 CLIFFORD S. PATLAK 


suitable analysis of experimental data and is not dependent on any 
particular transport model. 

The above conclusions may be derived from a consideration of 
equation (1) alone. This will show that equation (1) may be used 
as a general criterion for al] cases, not just the three examples 
mentioned above. If equation (1) is satisfied, then we see that 
certain flows must be linked together. For a molecule to be in 
equilibrium, its movement must not be totally dependent on the 
movement of any other molecule or molecules. ‘Therefore if the flow 
of a molecule is linked to the flow of other molecules, it cannot 
achieve a true equilibrium. (By linked, we of course do not mean 
such things as the effects of diffusion drag forces between mole- 
cules). 

We shall call all molecules whose coefficients in equation (1) 
are not zero linked-impermeative molecules. We may then define a 
permeative molecule as one which can traverse the membrane (izes 
is not impermeative) and which is not linked-impermeative. The 
last definition is the same as that given at the beginning of this 
section. 

From the above discussion we see that our attempt to decide 
whether a molecule is or is not actively transported by determining 
whether it is in equilibrium or pseudo-equilibrium would apply only 
to component molecules and not to linked-impermeative molecules. 
This is because in at least a pseudo-equiliorium state a linked- 
impermeative molecule will not in general be in equilibrium on the 
two sides of the membrane. Thus it would always appear to be 
actively transported on the basis of this criterion. Therefore, we 
limited our criterion to component molecules alone. 

It is, of course, possible that a linked-impermeative molecule 
will be actively transported. For example, in the more general 
exchange diffusion example discussed before, it is possible that 
@ will be actively transported from the internal to the external 
solution. Therefore R and/or S will appear to be actively transported 
from the external to the internal solution. However, equation (1) 
would still indicate that Q was a linked-impermeative molecule. 
For this case, as well as all other cases in which there is an active 
transport of a linked-impermeative molecule, we have no general 
criterion to indicate whether there is or is not active transport. 
The only thing to do is to thoroughly analyze the particular case. 
Then, if the mechanism leading to the linkage of the flows can be 


ACTIVE TRANSPORT 283 


determined, the decision of whether there is an active transport can 
be made by knowing the distribution of all of the molecules. Basi- 
cally this is because the distribution of the linked-impermeative 
molecules, even when the flows are zero, is determined both by 
the mechanism of the linkage and also by the initial distribution 
of all of the linked-impermeative molecules. 

The criterion of equation (1) may be exceedingly difficult to 
apply. We shall discuss three reasons for this. The first and major 
reason is due to the difficulty in making precise measurements in 
tracer experiments. However, since equation (1) may be written 
with all of the v’s as integers, then under the assumption of the 
‘simplicity of nature,’’ for most cases these integers will be 1, 2, 
or at most 3. Thus, if there is linked-impermeativity in the system, 
it may be detectable even though the values of the j’s may not be 
accurately determined. 

A second reason would be the difficulty, if not actual impossibility, 
in many experimental situations, of controlling both solutions of 
the system, as discussed previously in this section. In most cases 
this would involve the difficulty of determining the outflow of the 
different molecules and the difficulty of altering the concentration 
in the non-controllable solution and then knowing what the new 
concentration is. 

When the passive diffusion of an apparently linked-impermeative 
molecule is much smaller than the rate of crossing the membrane 
by means of the linking mechanism, we find a third case where the 
criterion of equation (1) might be difficult to apply. In this case, 
equation (1) could be satisfied within experimental error even though 
there really is no linked-impermeativity for this molecule. However, 
if this passive diffusion were small enough, then the system may 
appear to have all of its flows equal to zero at a time when sue 
molecule is not in equilibrium. Due to the finite time of the experi- 
ment, however, we may treat this situation as if this molecule 
were a linked-impermeative molecule and only investigate the other 
components. 

As a final comment in this section, we may note that it is possible 
for a permeative molecule to traverse the membrane by means of 
non-energetic carriers alone. As mentioned by T. Rosenberg (1948), 
this occurs when the molecule is carried across the membrane by 
two different carriers. It may also be accomplished by the method 
discussed earlier of allowing one of the group of molecules that is 
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carried by the more general exchange-diffusion carrier to also traverse 
the membrane by a different method. This also occurs if all the 
other molecules in the chemical compound method are permeative 
molecules. However, it may also be accomplished by means of a 
single carrier if, unlike the simple exchange-diffusion case, the 
carrier can itself traverse the membrane without the necessity of 
being bound to the molecule in question as in the example of the 
non-energetic carrier discussed previously. It is for this reason 
that we defined the term diffusion as we did at the beginning of 
this section. 

The remainder of this paper will be confined to systems of strong 
electrolytes in the pseudo-equilibrium state. 


Electroneutrality: 

It is commonly assumed in all works on the diffusion of electro- 
lytes that there is electric neutrality for all regions. The reason 
for this is perhaps best illustrated by E. Guggenheim (1950), when 
he shows that an excess of positive over negative ions as small as 
10” '° moles can lead to potentials of the order of 107 volts. In this 
section, potentials shall always signify electric potentials. The 
problem of electric neutrality in the study of diffusion of ions across 
a membrane has been considered and debated by M. Planck (1890), 
K. Sitte (1934), and J. Hermans (1936). It will be useful to recon- 
sider this problem, since certain rather interesting results can be 
found, 

For a potential difference to arise, there must be a net charge 
present somewhere. This may be explicitly expressed by the 
Poisson equation, 


2 4a 
VV= eae (2) 

in which V is the potential, p is the net charge density, and «x is 
the dielectric constant. For simplicity we shall assume that only 
univalent ions are present in the solution. (The general case may 
be treated in a completely analogous manner and similar conclusions 
will be reached), For an aqueous solution where C+ denotes the 
total positive ion concentration and C™ the total negative ion con- 
centration in moles per liter, the above equation reduces to 


2 
VV a 10 TCR eae (3) 
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For a membrane of at least 1u in thickness and a potential differ- 
ence across it of the order of 100 mV, we see that VV is of theorder 
of, or less than, 100 x 10-8 volts/1074 cm = 103 volts/cm. Since 
the potential is constant in the solutions on either side of the 
membrane and only changes in the region about and within the 
membrane, the gradient of the potential will change from zero to its 
average value in at most the width of the membrane. A better 
approximation would be that it has its maximal change in approxi- 
mately 1/10 of the membrane width or less (see discussion on 
J. Bartlett and R. Kromhout below). Therefore, for this case the 
total change of the gradient occurs within at least 0.1. Then 
Fin 10°/10-° volts/em = 10° volis/em. Thus, from equation 
(3) we have that (G— C) = 10-° mol/liter. Since biological 
concentrations may be > 10 ? mol/liter, we see that for the case of 
membranes 1y thick, and larger, the concentration difference between 
the positive and negative ions at a point is negligible compared with 
their actual concentrations. Putting it a different way, we see 
that for Poisson’s equation to be valid for ‘‘thick’? membranes, we 
may effectively assume electroneutrality, since a negligible deviation 
from this neutrality can account for any variation in the potential. 

If the membrane is ‘‘thin,’’ however, the conclusion of electric 
neutrality does not hold. For example, the membrane of a nerve 
fiber is usually taken to be of the order of 100 A. Using the same 
approach as above, we see that V?V =~ 10'? volts/cm?. Thus, 
(C* - C-) = 107! mol/liter which may actually be larger than the 
concentration of some ions and is of the order of the total concen- 
tration of the ions within and outside of the nerve fiber (A. Hodgkin, 
1951). 

If we had assumed that the membrane was composed of lipids, for 
example, and that the ions diffuse across such a membrane by pass- 
ing through the lipids instead of through aqueous pores, the ex- 
istence of lipids would lower the dielectric constant by a factor of 
about 30 within the membrane. However, the concentration of the 
ions in the membrane would also be very different from concen- 
trations found in the solutions about the membrane. This would be 
due to the very small partition coefficient of the lipids with respect 
to water. Thus the same conclusion as above would be reached. 

Let us now consider the situation first discussed by D. Goldman 
(1943), in which the gradient of the potential is constant across 
the membrane. Thus the gradient changes from a value of zero to 
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its maximum value in a distance of a few Angstrom units, and there- 
fore (C* — C-) would be of the order of 1 mol/liter. As pointed out 
by J. Harris (1956), this would cause a ‘‘condenser membrane’’ with 
charges on the surfaces, This effect may also be seen in the paper 
by J. Bartlett and R. Kromhout (1952), where Poisson’s equation 
was used to analyze the potential change across an infinitesimally 
thin membrane with a Donnan equilibrium in the surrounding solutions. 
In this case, they showed that the total change of potential effectively 
occurred within 100A; that the major change of the gradient of the 
potential occurred within 10 A; and also, of course, that the space 
charge was an important consideration. In any case, electric neu- 
trality is a questionable assumption when treating diffusion across 
a thin membrane (except, of course, under the Goldman hypothesis, 
and then only for the membrane proper). 

As seen above, we cannot rigorously assume electric neutrality 
at all points of the system for all cases. We shall, however, make 
this assumption for the remainder of this section because we will 
not deal with the membrane portion of the system. The assumption 
of electric neutrality is useful to qualitatively determine the presence 
and sign of a potential in a portion of a system. This fact is well 
known and has been implicitly used in the literature. However, we 
feel that this approach has not been exploited as much as it could 
be, and therefore we will discuss it at length. 

Let us consider a system in at least a pseudo-equilibrium state, 
which consists of a membrane separating two solutions denoted as 
the internal and the external solutions. The solutions may contain 
component ions, fixed charges, linked-impermeative ions, and imper- 
meative ions. We will assume that the partition coefficient is the 
same for both regions. Since we are in at least a pseudo-equilibrium 
state, all of the component ions are in equilibrium except for those 
ions which are in pseudo-equilibrium. Therefore in the external 
Solution the net charge of the component ions in equilibrium is 
equal and opposite to the net charge of the fixed charges, the 
linked-impermeative ions, the impermeative ions, and the ions in 
pseudo-equilibrium. If there were no potential differences between 
the internal and external solutions, then the concentration of the 
component ions in equilibrium would be the sare in the two solu- 
tions. (Actually, for equilibrium the activities, not the concentra- 
tions, are equal. However, since most biological solutions are di- 
lute and since the internal and external concentrations are approxi- 
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mately equal, the activity coefficients, calculated by means of the 
Debye-Huckel theory, are approximately equal in the two solutions. 
To very good approximation the concentrations may therefore be 
taken as equal). If the net charge of the fixed charges, the linked- 
impermeative ions, the impermeative ions, and the ions in pseudo- 
equilibrium, in the internal solution, is equal to the net charge of 
the same group of ions and charges in the external solution, elec- 
tric neutrality will then exist in the external solution. However, if 
the above net charges are noé equal, electric neutrality does not 
exist in the two solutions since the concentration of the compo- 
nent ions in equilibrium is the same in the two solutions. There 
must then be a potential difference between the internal and ex- 
ternal solutions for electric neutrality to exist. For example, if 
there were no fixed charges, linked-impermeative and impermeative 
ions, but one of the components is actively transported, then, since 
its concentration will be greater in the solutions toward which it is 
being transported, there will exist a potential difference between 
the two solutions. 

To find the sign of the potential difference we use the fact that a 
high potential region attracts negative ions from and repels posi- 
tive ions to a low potential region, and vice versa. Let the net 
charge of the fixed charges, the linked-impermeative ions, the 
impermeative ions and the pseudo-equilibrium ions in the in- 
ternal solution, be more positive than the net charge of the same 
group in the external region. Then in order to have electric neu- 
trality, more negative and fewer positive component ions which are 
not in pseudo-equilibrium are needed internally than externally. In 
order to accomplish this the internal potential must be positive 
relative to the external solution. On the other hand, if the internal 
net charge of the fixed charges, the linked-impermeative ions, the im- 
permeative ions, and the pseudo-equilibrium ions 18 more negative 
than the external net charge of the same group, then the internal 
potential must be negative relative to the external potential. er 
Thus, in a Donnan equilibrium, auiiete the xed charges oe ae 
impermeative ions are only in the internal coe the sign oO i 
internal potential relative to the external region is the same as e 
sign of the net charge of the fixed charges plus the impermeative lons. 


For the case of no fixed charges, linked-impermeative ions, and imper- 
the component ions is actively trans- 
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ported, the sign of the potential in 
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transported relative to the potential of the other region must be the 
same as the signofthe ion. Hence, we see that the potential difference 
resulting from the active transport of an ion will always oppose the 
active transport, as would be expected. There is also no difference 
in the sign or value of the potential difference between a fixed 
charge, an impermeative ion, or an actively transported ion. This 
can be seen from the derivation of the sign of the potentials and 
from the examples given above. Thus, we may say that pseudo- 
equilibrium is a type of Donnan equilibrium. 

The following is an example of the above approach applied to 
nerve fibers. If Na’ is the only ion actively transported outward 
across a nerve membrane and if its net charge outside is greater 
than the net charge of Na® plus the net charge of all the fixed and 
impermeative ions inside the nerve, then the outside potential is 
positive relative to the inside. This assumes, of course, that the 
nerve is in a pseudo-equilibrium state. As has been stated, this is 
not true for the nerve in vitro but may be true for the nerve in vivo. 

The following two factors relative to the above discussion should 
be made explicit. We will first discuss the consequences of a non- 
uniform charge distribution within a region. For example, it may be 
argued that most of the fixed charges in the internal solution occur 
at the boundary between the membrane and the internal solution. 
Since in pseudo-equilibrium j(X,2) = 0 everywhere and since we 
have electric neutrality everywhere, we may use the same argument 
as above for the potential at each point within a solution relative 
to the other solution. Thus, the above discussion regarding the Na* 
would apply to the potential in the axoplasm of the nerve away from 
the membrane, where the amount of fixed charge and impermeative 
ion would apply to the axoplasm and no¢ to the membrane. 

The second factor is concerned with what is meant by a fixed 
charge. If there is a bonding between one of the ions and the fixed 
material of the solution, then the resultant charge of the fixed 
material plus bound ion is part of the fixed charge. The concen- 
tration of the ion in solution will therefore be the total concentration 
of the ion (bound and free) minus the bound concentration. The 
concentration of the bound ion may be difficult to determine ex- 
perimentally. However, it can be readily calculated when the ion 
is known to be in equilibrium in the system; when the potential 
difference, the activity coefficient, and partition coefficient are 
known at all points; and when the bound ion concentration of one of 
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the solutioms is known. The same argument applies to the case of 
a bonding between one of the components and an impermeative 
molecule, where the combination is also impermeative. 


The general approach of using electroneutrality as a tool for 
demonstrating the existence and sign of a potential difference can 
be applied to the problem of active transport across a frog skin. To 
demonstrate whether or not an ion undergoes active transport, 
H. Ussing and K. Zerahn (1951) devised the clever approach of 
short circuiting the frog skin, This was accomplished by making 
the potential difference across the skin zero and making the concen- 
tration of all ions the same on both sides of the skin. This ex- 
perimental approach has also been used by L. Kirschner (1955). 
In explaining his results on the active transport of Na’, he assumed 
that the potential difference is zero at all points in the membrane. 
That this assumptionis actually false can be seen from the follow- 
ing argument. 

As shown by H. Ussing (1949a), active transport of Na* only 
occurs over a certain portion of the thickness of the membrane. 
The remainder of movement of the actively transported Na’ ion 
occurs in the rest of the membrane by passive diffusion. Consider 
a point of the region inside the membrane, between the inner solution 
and the active transport region, and immediately adjacent to the 
active transport region. The Na’ diffuses from this point to the 
inner solution. In the steady state, under this experimental condition, 
all other ions are in equilibrium. Then, if the potential difference 
between this point and the inner solution is zero, the diffusion flow 
of Nat must be due solely to a concentration gradient. But this 
implies that the concentration of the Na* ion is greater at this point 
than in the inner solution. Since the other ions are in equilibrium, 
their concentration is the same at this point as it is in the inner 
solution. By the same type of argument as before, there must be a 
potential difference between this point and the inner solution if 
electroneutrality is to hold. The same argument will, of course, 
hold for all points inside the membrane which do not lie within the 
active transport region. Whether or not there would be a potential 
difference present in the active transport portion would depend upon 
the particular mechanism of transport. 


Certain implicit assumptions were made in the above argument 
which will now be made explicit and discussed, One is that the 
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partition coefficients are the same at all points of the system. The 
above ‘argument will hold, even if this assumption were false, ex- 
cept for a highly improbable set of values for the partition coef- 
ficient. The fact that varying partition coefficients can cause a 
potential difference for any situation has been discussed in an 
analogous manner by M, Polissar (in F. Johnson, H. Eyring, and 
M. Polissar, 1954, p. 526). The same type of argument will apply 
when there are fixed charges in the membrane. 

Another point to be clarified is the assumption that the rate at 
which each Na” ion diffuses in the non-active-transport region is 
not much greater than its rate of active transport. If this were not 
true, the Nat inside the membrane would diffuse out at such a rapid 
rate that its concentration would be the same as the concentration 
of the internal solution, i.e., the effective boundary of the membrane 
would be at the boundary of the active transport region. In this 
case, although a potential difference is present, it is so small that 
for all practical purposes it may be ignored. The aoe result is 
true for the case in which the concentration of the Na’ is much 
smaller than the sum of the concentrations of the other ions. 


Types of Active Transport: 

In order to have a net active transport, the number of actively 
transported particles of any component transported in one direction 
per unit area and time must be greater than the number of particles 
of the same component actively transported in the opposite direction. 
Let A be a specific component which is actively transported. For 
the transport mechanism to be effective, A should be transported 
across a region of sufficiently low permeability so that the rate of 
passive diffusion is low. As pointed out by A. Krogh (1946), to 
bring about a pseudo-equilibrium state, passive diffusion must 
compensate for the active transport. Thus the greater the rate of 
active transport of A, and/or the smaller the permeability, the 
greater the difference in concentration on either side of this region 
must be to attain pseudo-equilibrium. As discussed in the previous 
Section, in pseudo-equilibrium the electric potential difference 
across the active transport region is a monotonic function of the 
concentration difference across the barrier. Thus the concentration 
difference alone and not the electric potential difference is all 
that is really necessary to consider in any qualitative discussion 
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of the effect of different active transport rates, potential barrier 
heights, ionic mobilities, partition coefficients, etc. 

Consider a system in pseudo-equilibrium. Let the active transport 
of all A’s take place across a region which is not the entire membrane, 
and assume that there are no linked-impermeative molecules in the 
system. Then all of the non-actively-transported components will 
be in equilibrium at all points of the membrane. The 4’s will also 
be in equilibrium at all points on either side of the active transport 
region, but, of course, not across this region. One of the major 
characteristics of equilibrium is that the concentration at one point 
relative to another point depends only upon the characteristics of 
the two points. From this we see that the characteristics of non- 
actively-transporting regions of the membrane do not play any role 
in determining the pseudo-equilibrium concentrations and the electric 
potential distributions in the two solutions, except insofar as the 
impermeativity and linked-impermeativity characteristics of certain 
molecules may be due to these regions. They only determine the 
rate at which tracers flow across the membrane and the rate at which 
the system approaches equilibrium when the concentrations of the 
components in the internal and/or the external solutions are modified. 

As stated in the Introduction, if one of the two solutions of a 
system cannot be controlled, it becomes difficult to differentiate 
between the true and pseudo-equilibrium states. By controlling a 
solution, we of course mean that all of the properties and parameters 
of the solution are known and can be altered at will. Examples of 
non-controllable solutions can be found in most biological systems 
in which an inner solution is completely surrounded by a membrane, 
such as muscle fiber, mitochondria, and erythrocytes. For such 
systems a method which could be used to determine whether the 
system is in true or pseudo-equilibrium is the application of a drug 
to the membrane which alters the properties of only the membrane 
(without allowing the linked-impermeative and impermeative mole- 
cules to become permeative). This is done either by changing the 
membrane permeability and/or affecting the active transport mecha- 
nism. After application of this drug to the membrane let the syetom 
attain a new true or pseudo-equilibrium state. Changes in the 
concentration distribution of any or all of the components in the two 
solutions would then indicate that the original system was in the 


pseudo-equilibrium state. This is based on the fact, discussed 


292 CLIFFORD S. PATLAK 


above, that if the original system were in the equilibrium state, its 
component distribution is independent of the membrane properties. 
Therefore, altering the membrane should not alter the component 
distribution. Conversely, if the concentration distribution does not 
change, this would indicate that the system was either initially in 
the equilibrium state or that the drug did not affect the active trans- 
port region, although it may have affected the membrane in the non- 
active-transport region. It is absolutely essential, however, that 
this drug does not affect the internal or external solutions but only 
affects the membrane. This is a stringent requirement because of 
the difficulty in determining whether a drug affects the two solutions. 

Active transport may be accomplished in any one or a combination 
of two general ways. One is via a mechanism which acts on all of 
the A’s simultaneously and shall be here designated a general type 
of active transport. This may be accomplished by either a push 
(as in a water flow carrying ions with it) or by a pull (as in a dif- 
fusion drag force, such as is found in the examples in the Intro- 
duction). 

The second active transport mechanism is one which acts on 
one or a few molecules of A at a time and which we shall call an 
individual type of active transport. This may be accomplished by 
either carrying A across the active transport region (by means of a 
carrier) or by allowing the particle to go across by ttself under 
simpler conditions. Thelattermay be accomplished by: (a) decreasing 
the resistance of the region (by decreasing the activation energy of 
the region to flow, or by opening a “‘hole’® and thus allowing A to 
take an effective ‘short-cut’? across the region); (b) increasing 
the energy of A (by a chemical reaction which then allows A to 
pass across the active transport region more readily and, also, 
perhaps more quickly if the energy increase causes a translational 
velocity increase); or (c) altering the properties of A (by forming 
an isomer of itself which can pass across the region more readily). 
The mechanism allowing the particle to go across by itself under 
the above simpler conditions will be mediated by a ‘*Maxwell 
Demon.’ This will be discussed at greater length in the next 
section. It is also possible to have a negative Maxwell Demon 
which prevents A from crossing the membrane and therefore acts 
somewhat like an active transport in the opposite direction. 

Some nomenclature and equations to be used in the remainder of 
this paper will be introduced before discussing the two types of 
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active transport. We shall let C denote the concentration of the 
component. The first subscript, which will be & or p, denotes the 
Ath component which is in equilibrium or the pth component in 
pseudo-equil’brium, respectively; the second subscript, which will 
be 2 or e, denotes the internal or external solution respectively; 
the superscript, which will be + or —, denotes a positive or negative 
ion; Z denotes the absolute value of the valence of the ion; V 
denotes the value of the electric potential at a point relative to the 
internal solution; w denotes the mobility; R, 7, and F denote the 
gas constant, the absolute temperature, and the Faraday constant, 
respectively; denotes the electrochemical potential; 1, denotes 
the standard chemical potential; 6(«) = exp{-[y, (2) - 4, ,J/RT} 
denotes the partition coefficient of the point 2 relative to the inner 
solution (F. MacDougall, 1939); 7 denotes the net flow of a component 
across a unit area per unit time; and y denotes the activity coef- 
ficient. 

If we assume that the molecule is flowing due to pure diffusion 
and that there are no diffusion drag forces between the molecules, 
then (L. Onsager and J. Fuoss, 1932) 


+ wea eee 
j=- uC, Va. (4) 


When the above equation is expressed in rectangular, cylindrical, or 
spherical coordinates, and the concentration is a function of only 
one of the rectangular coordinates or only the radius of the cylindrical 
or spherical coordinates, then all three coordinate systems give the 
same equation. In the following, z will denote the coordinate along 
which the concentration varies. Then equation (4) reduces to 


+ on 
k ke ek da 
As is well known, 
ray pe +RT log Ci +RT log ye Lega. (6) 


Inserting the above into equation (5), we have, after slight rear- 


rangement, 
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(7) 


The first term on the right, involving the dCt/de, has been dis- 
cussed extensively by L. Longsworth (1933) for a slightly more 
limited case. For equilibrium, the above reduces to a modified 
Maxwell-Boltzmann distribution, 


+ 
bee 
ye 4 ky 
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In the literature, yy is usually assumed equal to unity (although all 
that is really needed for the equilibrium case is that the y’s be 
the same at all points). Also, 6, . iS assumed equal to unity (i.e., 
the solvent is the same in the two solutions considered). As is 
well known, equation (8) is a thermodynamic result. It is therefore 
valid even though equation (4) may not be applicable (cf. the 
passive diffusion of K+ through nerve fiber, A. Fodgkin and R. 
Keynes, 1955). 


1. General type of active trans port: 

The mathematics of this subject has been extensively treated 
in the literature (cf. H. Linderholm, 1952, for perhaps the most 
complete treatment). Therefore we will not deal with it at length. 
Basically, the general type of active transport implies an ‘“‘active- 
transport’? force applied to the molecules in question. Letting 
SEZ Ob3/ dx represent this force (H. Landahl, 1953; H. Ussing 
and Kk. Zerahn, loc. cit.) in which by = 0, the basic equation for 
this type of active transport is 


on * dbp 
ot tpt 
j-upot| Me a pzs be e (9) 


For equilibrium, the above equation may be solved for oh - to 
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yield 
Fz* ‘ 
+ a = Veto 2) 
® iad Pp,e 
P,e 


Since Ger Ct yeh: ye ono and V, may be evaluated ex- 
perimentally, the active transport term (b> -) may be determined. 
We may note that equation (10) contains no explicit terms involv- 
ing the membrane. 

We are familiar with only two basic models for the general type 
of active transport, although there are probably many others. One 
model is a water flow which drags A along with the flow (H. 
Ussing, 1949b). The water flow usually arises from a difference in 
osmotic pressure across the membrane, although it may arise from 
a@ pumping action (as in the circulatory system). The major dif- 
ficulty with this model is that the water flow would drag all of the 
permeative ions with it and would not be specific for any partic- 
ular ion. 

The other model arises from the existence of diffusion drag 
forces. In this case a passively diffusing material drags the 
actively transported substance. This drag is accomplished by 
some unspecified mechanism, such as an attractive force of some 
type, or actually colliding and pushing the actively transported 
molecule (cf. N. Rashevsky, 1948). (It should be emphasized that 
the ‘“‘carrier’? model is not a satisfactory model for this type of 
transport, as discussed below). Note that both the above models 
possess the common feature that the material which is doing the 
‘forcing’? (pushing or pulling) cannot be in a true or pseudo- 
equilibrium: state but must be either in a steady or a transient 


state. 


2. Individual type of active transport: 


Although almost all models which have been invented to ex- 
plain the different experimental cases of active transport are of 
this type, there have been no mathematical treatments of it in its 
most generalized form. Therefore, it might be of interest to ex- 
plore the mathematics of this situation. 

When considering the general type of active transport, we as- 
sumed a force on all of the actively transported molecules. This 
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force was superimposed upon all other forces acting on the mole- 
cules (such as the electric potential gradient, etc.), and this force 
caused a net ‘‘drift?? of each molecule. However, a molecule is 
free to move back and forth, as it does under pure diffusion, but 
due to the action of this force, it would move statistically more in 
one direction that in the other. Thus we had a net movement of 
the particle due to the diffusion force, the force of the electric 
potential, and this active transport force. All three forces were in 
series. Therefore the net flow was the sum of the flows due to 
each of the separate forces alone. Since the active transport 
force independently acts on each particle, we saw that the net; 
flow due to this force alone was Fux Cx F Zs ) o3/ dx. 

The situation for the case of the individual type of active trans- 
port is different. (The following discussion may be clearer if the 
reader keeps in mind the carrier model as a particular example of 
this type of active transport.) Here, A is acted upon, or “picked 
up,’’ at a certain point and transported by some means to a second 
point. At this point it is allowed to re-enter, or is ‘‘discharged,’” 
as @ passive (i.¢., non-actively-transported) molecule. When the 
molecule is being transported it is usually not allowed to free 
itself from the transporting mechanism and become a passive mole- 
cule. Also, passive molecules in the region between the points 
of “‘pick-up’’ and ‘‘discharge’? are usually not allowed to be trans- 
ported by the active transport mechanism. [This non-interaction 
is introduced because almost all models with which we are famile 
iar have this property, although it is, of course, possible that 
there will be the type of interaction mentioned above. An example 
of this type has been treated by L. Kirschner (loc. cit.). However, 
in the mathematical derivations to follow, we will assume that 
there is no interaction between the passive and actively trans- 
ported A.} Thus, the flow of A due to the active transport will 
not be proportional to the concentration of A at each point, as it 
was in the general type of active transport, but will be a constant 
value determined solely by the properties of the membrane, the 
electric potential distribution, and the concentrations at the points 
where A is first acted upon and where it is released. Hence the 
flow due to this mechanism can not be represented by a force term, 
as in equation (9). Also, the active transport flow is the same at 
all points within the region of active transport and is, of course, 


ACTIVE TRANSPORT 297 


zero at all points outside this region. Since under the assumption 
that there is no interaction between the actively transported A’s 
and passive molecules, the passive molecules in the active trans- 
port region will diffuse as if there were no active transport (i.e., 
active transport and diffusion flow are in parallel). 

If we let J, represent the net flow across a unit area per unit 
time due to the active transport mechanism, the basic equation 
becomes 


(11) 


in which J, is constant within the active transport region and zero 
outside of it. Substituting into the above equation the value for 
Th from equation (6), we have 


oct ) : ov 
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In equilibrium 1s = (0. Equation (11) then reduces to 
(13) 


in which the mobility, us cannot be cancelled out. Thus, the con- 
centration at all points will be a function of the parameters of the 
system (in particular those of the membrane), even though the 
active transport system (represented by J5) may be independent of 
the membrane properties. This is contrary to the results found in 
the general type of active transport. 3 

To illustrate some properties of the pseudo-equilibrium state we 
shall solve the above equations for a simple case and then discuss 
the assumptions involved. This will provide additional insight 
into some of the problems connected with this approach. We shall 
assume (a) that there is only one actively transported ion (which 
we denote by the subscript 7); (b) the valence of A is + 1; (c) the 
boundaries of the active transport region (i.¢., that portion of the 
membrane where J, 4 0, and which will be denoted by the second 
subscript a) are at 2 = and 2 = B; (d) all of the passively diffus- 
ing ions are univalent (Z 5 = 1); (e) ifere are no linked Nppermealive 
ions; (f) J, is positive (ie., the active transport 1s from the in- 


298 CLIFFORD S. PATLAK 


ternal to the external solution); (g) there is electroneutrality at all 
points; (h) there are no fixed charges within the active transport 
region; (i) is Yp Dry bos and wu, are constant within the active 
transport region; (j) the net charge of the fixed charges and imper- 
meative ions in the solutions will be denoted by B; (k) while being 
actively transported, A plus the transporting mechanism, as well 
as all of the non-component molecules in the membrane, are elec- 
trically neutral; and (2) the internal and external solutions are 
uniform in character. 
In pseudo-equilibrium equation (7) then becomes 


DOR 0 ST ike OV 
26s cos at fog? ed 
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and equation (12) becomes 


dC ) FE .0¥ J 
Fi 0,| Poet + zoe, (15) 
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From equations (14) and (15) we now add all of the equations for 
positive ions and subtract from their sum the sum of the equations 
for negative ions. Using the assumptions of electroneutrality at 
all points, and the constancy of the y’s and 6’s within the active 
transport region, we obtain for the region of active transport 


2F oV J 
—— —— Cy = pP__, 
RT de 2  Ractiak ted Ses 
Solving equation (14) we obtain for the active transport region 
# a8 
VE + — V(x) 
Ci (@) oy Cpe See oaks (17) 
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Inserting equation (17) into equation (16), and simplifying, we obtain 
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In many experimental cases, es can be easily measured. For the 
case in which it is Ce instead, which can be most readily meas- 
ured, the results derived below can be transformed quite easily 
into the desired form. Thus, we shall attempt to solve for V, in 
terms of the parameters of the membrane and Coe Equation (14) 
may be solved to yield 


+ + - _F = 
Vi (2) bre SF ep Me VO 


C7. =C, (2) 
k,e k P 
Vic , @ 


(19) 


Since there is no flow or active transport across the boundaries be- 
tween the active transport region and the non-active-transport 
regions, the A inside of and at one boundary of the active transport 
region is in equilibrium with the A in the non-active-transport 
region on the other side of that boundary. A similar equilibrium 
exists across the other boundary of the active transport region. 
Equation (15) may then be solved for the regions outside of the 
active transport region to yield 


- — [y,-vA)] 
Yp,a One RT A 
Cie = ©,(B) aa é ’ (20) 
Pre p,a 
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Assuming electroneutrality and using equations (19) (with 2 = B) 
and (20), we can solve for exp {FV(B)/RT} in terms of the con- 
centration in the external solution, V,, and the partition and ac- 
tivity coefficients of the active transport region and the external 
solution. Similarly, by use of electroneutrality and equations (17) 
(with 2 = «) and (21), we can solve for exp {FV (a)/RT} in terms 
of the concentrations in the internal solution, and of the partition 
and activity coefficients of the active transport region and the 
internal solution. Then by use of electroneutrality and equation 
(8) we have, after lengthy but elementary manipulation, the fol- 


lowing implicit equation for V,: 
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The above equation is too complicated for practical use. There- 
fore, the following simplifications are made. Assume (m) that the 
net charge of the fixed charges and impermeative ions in both the 
internal and external solutions is zero; (n) that the activity co- 
efficients at each point are the same for all of the different compo- 
nents; and (0) similarly for the partition coefficients. That is; 
Vie (2) = yp (x) = y(a) and bf (x) = b, (2) = b(z). The above equation 
then reduces to 


2 RT 6 Pe j a 
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Since Uy ag 0, %T/y,(B - %) is the permeability of the active trans- 
port region to A (which we denote as Pig) (cf. A. Hodgkin and 
B. Katz, 1949), equation (23) reduces to 
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The value of V, in terms of the internal region concentrations may 
be obtained under the same assumptions used in deriving the above 
equation by inserting equation (8) into equation (24) and rearrang- 
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ing. We then obtain 
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Therefore, to increase V, we may either increase Jy, decrease 
P,,a, Or decrease the total ion concentration. This would be 
expected on the basis of the discussion at the beginning of this 
section and the discussion on the sign of the potential for the 
transport of a positive ion. Eowever, by changing the concentra- 
tion of ions and measuring the resulting V, we can only find the 
ratio J,/P, q, and neither J, or P,, individually. Contrary to the 
results found for the general type of active transport [equation 
(10)], the resultant electric potential expressed in the above equa- 
tions does depend on the properties of the membrane. This agrees 
with the conclusions made from equation (13). Specifically, the 
resultant electric potential depends only on the properties of the 
active transport region and does not depend on the properties of 
any other portion of the membrane, in agreement with the state- 
ments at the beginning of this section. Although it is probable 
that P,, is constant as the concentration of the ions is changed, 
J» probably is not. In fact, J, is most likely a function of C (4) 
and V(z) and possibly of many other variables. Thus, by changing 
the concentration of the ions, which will cause a change of both 
C ,(%) and V(z), J, will change in a manner which cannot be deter- 
mined from equation (24) alone, for not only is P, , not known, 
but the values of C ,() and V(z) are not usually measurable. 

We will now discuss the assumptions involved in the above 
derivation in the order in which they were presented, beginning on 
page 297. 

The assumption (a) of only one actively transported ion appears 
to be valid for many actual cases. However, there are many 
cases of “‘linked’? transport in which two components are trans- 
ported in a stoichiometric ratio. Also, there are many cases in 
which a number of components are actively transported independ- 
ently of one another [e.g., in marine alga (G. Scott and H. Haywood, 
1955) and in erythrocytes (M. Maizels, 1954)]. The two cases 


302 CLIFFORD S. PATLAK 


above can be treated in a manner similar to the treatment with as- 
sumption (a). 

The assumptions (b) and (d) of univalent ions and the assump- 
tion (e) of no linked-impermeative ions are introduced solely for 
mathematical simplicity. If they were not used, the problem be- 
comes extremely difficult to treat. 

The assumption (f) that J, is positive merely specifies the 
internal and external solutions. That is, for nerve the external 
and internal solutions would be the outside solution and the inner 
core of the nerve fiber respectively. However, for frog skin, the 
external and internal solutions would correspond to the inside and 
outside respectively. 

The assumption (g) of electroneutrality may be considered valid 
for a wide active transport region. From the discussion of the pre- 
vious section, however, if the active transport region is thin, as 
appears to be the case for peripheral nerve, erythrocytes, and frog 
skin, the assumption of electroneutrality may not be valid. This is 
expecially true since the partition coefficient for the active trans- 
port region may be expected to be very small. Hence, in this case 
if the equations are to be solved exactly they must be solved with 
the use of the Poisson equation, which introduces a considerable 
mathematical complexity, 

The assumption (h) that there is no fixed charge in the active 
transport region was used for simplicity. The case of a fixed 
charge in this region could easily be treated in a manner analogous 
to the discussion above. 

The assumption (i) of constant activity and partition coefficients 
and the non-varying mobility of A within the active transport region 
is probably valid for most biological situations in which the active 
transport region is very thin, However, hidden in this assumption 
is a rather restrictive condition. If we consider a specific homo- 
geneous physical region in which the active transport occurs and 
if A is *‘picked up’? and ‘“‘discharged’® within this region, then the 
above assumption will be valid. However, A may be ‘‘picked up’’ 
from or ‘‘discharged’? into the regions immediately outside of the 
above defined physical region. (This may perhaps best be illus- 
trated by the example of the entire membrane being the active 
transport region and the carrier attaching itself to A in the internal 
solution rather than in the membrane itself), In this case, the 
active transport region would have to include that portion of the 
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system immediately outside of the physical region in which 4 is 
either ‘‘picked up’? and/or ‘‘discharged’? (since J, # 0 across the 
boundary of the physical region). Thus, it would then be physi- 
cally improbable that the assumption of constant parameters would 
be valid. Further discussion of this assumption will be found 
later in this section. 

The assumption (j) of the presence of fixed charges and imper- 
meative ions in the internal solution was introduced for the sake of 
generality. However, in order to avoid a quadratic equation they 
must be absent, as in assumption (m). 

The assumption (k) of the neutrality of all non-component mole- 
cules in the membrane, as well as the neutrality of A plus the 
transporting mechanism while A is actively transported, is obvi- 
ously false. If A were actively transported without being neutral- 
ized, then, since it is charged it is not neutral. If A were neutral- 
ized while being actively transported, there must be charged non- 
component particles in the system. (This must be so since elec- 
tric charge cannot be created or destroyed, since there is only one 
actively transported ion, and since, for the entire system to be in 
the pseudo-equilibrium state, the particles in the active transport 
region must be in at least the steady state.) This may be made 
clear if we consider the example of the diffusible carrier system. 
If A were neutralized by the charge on the carrier, then the carrier, 
upon diffusing back alone, would be charged. 

In certain cases however the assumption of neutrality is a good 
approximation. We shall mention two general examples. The first 
is that in which active transport takes place through a portion of 
the membrane which is different from the region of passive trans- 
port. This may be accomplished by having the active transport 
take place in special channels separate from the channels used 
for passive diffusion. (This mechanism will be discussed ex- 
plicitly later in this section.) In this case, the equations which 
have been used are still applicable to the passive diffusion region. 
Thus, if the active transport channels are sufficiently isolated from 
the passive diffusion regions, equation (24) will be valid. The 
second example is found when the rate of movement of each mole- 
cule of A in the active transport state is very high compared to its 
rate in the passive diffusion state (all other particles, except for 
A in the active transport state, being neutral). (For the case of 
the carrier model, this would mean that the carrier-ion complex 
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has a mobility which is much greater than the mobility of A alone). 
The concentration of the passively diffusing A must be much greater 
than the concentration of the actively transported A. This is true 
since for pseudo-equilibrium the net rate of the passive diffusion 
of A must be equal and opposite to the net rate of the flow of A 
in the active transport state. Hence, with good approximation we 
may ignore the contribution of A in the active transport state to 
the equations of electroneutrality in the active transport region. 
(The presence of non-actively transported ions makes this argument 
even stronger.) Thus, the assumption of the neutrality of all other 
particles is a good approximation for this example. For a situation 
in which the charge on the actively transported A or on any other 
particle must be considered, a specific model of the active trans- 
port mechanism would have to be postulated. 

The assumption (/)-of the uniformity of internal and external 
solution is valid in many experimental situations. In reference to 
this, note the comment of the previous section concerning the dis- 
cussion of non-uniformity within a region. 

The simplifying assumption (n), page 300, that the activity 
coefficient is the same for all ions is valid for most actual cases, 
This is because the concentrations usually used in biological ex- 
periments are low so that the activity coefficients are approxi 
mately one. The activity coefficient of a low concentration of ions 
may be given by the Debye-Huckel theory. This theory gives 
log Vi = ~constant-(Zt)?, where the constant, which depends pri- 
marily on the ionic strength, is the same for all ions. A similar 
result applies to the Yp Since all of the ions are assumed to be 
univalent, all of the y’s are the same at any point. This also 
points up the fact that since the total ionic concentration differs 
at different points, the above constant differs at different points, 
Therefore Yi and Yp are not constant as a function of z, as we have 
assumed. However, since the concentrations are low and since 
the total ionic concentration may not change by much across the 
active transport region, the assumption of the constancy may be 
a good approximation. Also, due to the complexity of the re- 
lationship between y* and yp and the total ionic concentration, 
any attempt to give an exact treatment will immediately lead to 
insuperable mathematical difficulties, 

The simplifying assumption (0) that the partition coefficient is 
the same for all ions does not hold too well; although on the basis 
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of the solubilities of KCl and NaCl in various alcohols, it may be 
valid to within an order of magnitude. In any case, this assump- 
tion was necessary to discuss equation (22). 

There are two implicit assumptions worthy of discussion. The 
first is that the distribution of A is given by equation (15). In 
most cases equation (15) will be valid. One case in which it will 
not hold is that in which the channel through which the ion diffuses 
is narrow enough so that steric effects enter into the diffusion con- 
siderations (cf. H. Davson and J. Danielli, 1948; and A. Hodgkins 
and R. Keynes, loc. cit.). (The distribution of ions in equilibrium 
is always given by equation (14), since this equation is merely a 
restatement of the constancy of the electrochemical potential 
which is always true in equilibrium.) The second implicit assump- 
tion is that the diffusion can be treated by a one-dimensional 
equation. That this is not always correct is most evident in the 
case of the erythrocyte, for example, where it has been shown that 
only 0.1% of the surface of the erythrocyte is available for dif- 
fusion (A. K. Parpart and J. F. Hoffman, 1954). 

The assumptions have been discussed at great length to indicate 
that even in the simplest case that can be considered (i.e., pseudo- 
equilibrium which is much simpler than the case of steady state), 
there are a great many difficulties which confront us. We have also 
seen that most of the assumptions will not be valid when applied 
to many actual biological situations. However, equation (24) is 
probably a good approximation for many actual cases. 

We will now discuss in detail a model of separate membrane 
channels for active and passive diffusion mentioned earlier. In 
this model we will use the same simplifying assumptions as were 
used to derive equation (24). For this model, we further assume 
that all of the channels under each of the two different categories 
are identical, that there is no interaction between channels of the 
different categories, that both types of channels are of the same 
length, and that there is complete mixing in the regions immedi- 
ately adjacent to the boundaries of these channels. The latter 
assumption of complete mixing implies that the mobility of A in 
the regions adjacent to the two types of channels is much larger 
than in either of the channels. Of course, the active transport 
channel will be composed of a section in which there is active 
transport and, on either side of it, sections in which there is only 
passive diffusion. Let the subscript d refer to the passive dif- 
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fusion channel; the subscript ¢ to the active transport portion of 
the active transport channel; the subscript dt to the passive dif- 
fusion portion of the active transport channel, Let J,, be the flow 
per unit area and time through the active transport channel; J, be 
the flow per unit area and time due to the active transport mecha- 
nism, which is the same definition used before; 6 be the ratio of 
the total cross sectional area of all of the active transport chan- 
nels to the total cross sectional area of all of the passive dif- 
fusing channels. Let ¢ and 7 refer to the inner and outer boundary, 
respectively, of the channels, while & and f refer to the inner and 
outer boundaries, respectively, of the active transport portion of 
the active transport channel. Since we have assumed a pseudo- 
equilibrium state, the total transport of A through the active trans- 
port channels must be equal and opposite to the total transport of 
A through the passive diffusion channels. 

Therefore 
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Adding the last two lines above and subtracting from their sum the 
first line above, we have, after elementary rearrangement, 


J 
Jy = E 5 (28) 


[soo S i 5 
1 + ee + 
Boyd Pp,d 


We will further assume that the mobility of A across the boundaries, 
€ and 7, between the passive diffusion channel and the regions 
immediately adjacent to this channel, is much larger than the 
mobility, Una? of A inside the passive diffusion channel. Thus, A 
inside of and at one boundary of the passive diffusion channel will 
effectively be in equilibrium with the A in the region on the other 


side of this boundary. A similar situation will exist at the other 
boundary of the passive diffusion channel. We then insert equation 
(28) into the first line of equation (27), use the same method as 
has been used before, and make the following substitution 


P PaasP 
eee a (29) 


p,dt 


We now find our result to be exactly the same as equation (24). 
This indicates that when the above assumptions are made and the 
system is in pseudo-equilibrium, we cannot distinguish between 
the existence or non-existence of separate channels by measure- 
ments of only the internal and external parameters. 

To complete our results on the pseudo-equilibrium state let us 
calculate the concentration of A in the external region for a given 
V_, the total external negative ion concentration, and the internal 
concentration of A. We will assume that the system is in the 
pseudo-equilibrium state, that all ions are univalent, that the 
valence of A is +1, and that the activity and partition coefficients 
are the same for the internal and external regions. By use of equa- 
tion (8), and the assumption of electroneutrality, we have 
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Note that the above equation is valid irrespective of the particular 


mechanism of active transport. 
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As discussed at the beginning of this section, models for the 
individual type of active transport fall under two main headings: 
carrier-type models and Maxwell Demon models. We shall first 
consider the carrieretype model. In this model A becomes attached 
to the carrier, forming a complex of some type. By means of dif- 
fusion, rotation, or structural change of the complex, A is then 
transferred across a region (usually containing a potential bar- 
rier). On the other side of the region the complex breaks up into 
A and the carrier. Energy is supplied to the system in the reac- 
tion of the formation of the carrier-molecule complex and/or in the 
reaction of the break-up. Most of the theoretical work and model 
construction for active transport mechanisms has been concerned 
with this model. Hence we need only very briefly outline some of 
the general types used. 

There are three basic carrier models. The first model we may 
designate the fixved-carrier model. (This has been called the 
“‘reaction-type’? model by T. Rosenberg and W. Wilbrandt (loc. 
cit.). However, this term includes other models besides the fixed- 
carrier models.) In this model A attaches itself to the carrier. 
Without any translational movement of the carrier, the complex 
then transports A across a thin but high potential barrier and then 
releases A from the carrier molecule. The transport of A can be 
accomplished either by a rotation of the complex or by having A 
attach itself to the ‘‘free’? end of a long carrier molecule which 
then contracts if initially expanded or expands if initially con- 
tracted. 

The second model is a generalization of the previous one. That 
is, there are many carriers such that A joins onto the first car- 
rier and is transferred to the second carrier which then transfers 
it to the next carrier, etc. In this manner A is transported across 
the region (cf. H. Lundegardh, 1954). 

The third type of carrier model is the most extensively treated 
carrier model in the literature. In this model A combines with the 
carrier to form a complex which diffuses across the region and then 
Separates on the other side. In this case, the carrier may be a 
complex molecule; it may be a simple electron or proton (H*) which 
may act to neutralize an ion so that it can more easily cross a 
lipid layer; or it may be nothing as when the ion is transformed 
into an isomer of itself which can then diffuse across the region. 
The last case is included under the heading of a carrier mechan- 
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ism because it follows the same pattern as the other carrier wodels 
and is just a limiting case. However, it can also be included under 
the Maxwell Demon model discussed below. Examples of the 
third type of carrier model are extremely numerous and may be 
found in, or from, the other references listed in this paper. 

Without necessarily considering the effect of the electric po- 
tential, general theoretical treatments for the carrier model may 
be found in J. Franck and J. Mayer (1947), J. Best (1949), M. 
Polissar (loc. cit.), and T. Rosenberg and W. Wilbrandt (loc. cit.). 

We will now consider the Maxwell Demon model for the indi- 
vidual type of active transport. To the best of our knowledge it 
has not been extensively considered in the literature. However, 
there have been four allusions to the possibility of its existence. 
We shall discuss this model in the next section. 


Maxwell Demon: 

The Maxwell Demon is a hypothetical mechanism which allows 
each particle of A to cross a potential barrier by itself more read- 
ily in one direction than in the other. It differs from a carrier in 
that no part of it in contact with A moves while A is crossing the 
barrier. The Demon operates either by lowering the potential 
barrier height for A, by increasing A’s speed of crossing the bars 
rier,sor by modifying A so that it can cross the potential barrier 
more readily. 

The Maxwell Demon was first proposed by J. C. Maxwell (1916, 
p. 338) in order to suggest a method of violating the second law of 
thermodynamics. His major assumption was that the Demon will 
react specifically to one type of molecule and allow it to be trans- 
ferred more readily in one direction than in another across a poe 
tential barrier. So that the second law of thermodynamics will not 
be violated, energy must be supplied to the system and must be an 
integral part of any proposed model. This has been implicitly 
stated by R. Dean (1941), who recognized that a Maxwell Demon 
could account for active transport, but who then dismissed the 
idea on the basis of the need of an energy source. , 

For a Maxwell Demon to affect the movement of A it must inter- 
act with A. Since there does not seem to be any evidence for long 
range forces, it therefore must come into virtually direct contact 
with A or come in contact with some other particle that contains 
A. When the Demon reacts with A or with a complex containing A, 
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energy must be liberated in order to transfer A across a barrier. | 
This energy may either come from the compound which formed a 
complex with A or it may arise from a reaction of a substrate with 
the Maxwell Demon. In either case, after all of the reactions the 
Demon must be restored to its original state. While the Maxwell 
Demon is active and is allowing A to cross a potential barrier more 
readily, it may also (depending upon the particular type of Demon) 
allow a particle to cross in the opposite direction. However, there 
will be a net transport in the required direction for all cases. This 
is true since A always crosses in the specified direction, and the 
particle which is allowed to cross in the other direction may not be 
present at the required position at the necessary time, or the 
particle may not be an A. 

We shall primarily concern ourselves with Maxwell Demons which 
transfer A across a potential barrier. However, there may also be 
negative Maxwell Demons which may keep A from crossing the 
membrane. That is, a net active transport of A from the internal to 
the external region may be accomplished either by aiding the trans- 
port from the internal to the external region, hindering the trans- 
port from the external to the internal region, or both. If it is to 
hinder efficiently, the negative demon should be located at a posi- 
tion immediately beyond the highest potential barrier so that it 
need only react with a minimal number of A. For example, were it 
placed by the external solution it would react with many mole- 
cules of A which would not have gone through the membrane in any 
case. 

As discussed earlier, one possible mode of operation of a Max- 
well Demon would be for it to increase the velocity of A through a 
chemical reaction. This would not only increase the velocity of A 
across the potential barrier, but since its kinetic energy would be 
increased, it would allow it to pass over the barrier more easily, 
This idea has been suggested by O. Schmitt (1951, p. 86) and dis- 
cussed more fully by M. Polissar (loc. cit., p. 760). 

Besides the comments made here, and by the preceding two 
authors, on this model, there are two other factors which should be 
mentioned, First, there is the fact that A in solution is surrounded 
by its water of hydration as well as the water of the solution 
about it, and this will help determine a maximum barrier width. 
In general, if a particle is given an increase of velocity, it will 
lose this excess energy due to collisions with the surroundings in 
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_ the space of about 1A because of the “‘cage effect’? (J. Franck and 
E. Rabinowitsch, 1934, p. 125; W. Libby, 1947). Thus, if the 
mechanism of velocity increase due to a chemical reaction is to 
work, the reaction must occur at the boundary of the major barrier 
to be crossed, which must then be of the order of 1A in width. 
Eowever, the Maxwell Demon could be placed inside the membrane, 
where the distance between collisions might be slightly larger than 
1A. However, it is hard to imagine that the figure of 1A could be 
increased by very much. 

The second factor concerns the efficiency and effectiveness of 
the active transport mechanism. If the system is to be effective, 
the permeability of the barrier, P, ,,:must be small. This in- 
dicates that the activation energy of crossing the barrier, AE i 
must be large, since P, , is proportional to exp {+ AE '/RT} (S. 
Glasstone, K. Laidler, and H. Eyring, 1941). If the active trans- 
port and the passive diffusion occur by means of separate chan- 
nels, then the above statement is even stronger, as can be seen 
from equation (29), If AE lis of the order of or greater than the 
energy given to A by the reaction, A will still have a certain 
probability of passing over the energy barrier, due to the Maxwell- 
Boltzmann distribution of energies of A when it is bound to the 
Maxwell Demon. This probability is, of course, much larger than 
if there were no reaction, but it may not be high enough to ac- 
count for the efficiencies that are experimentally observed. 

Another model, which was independently proposed by Professor 
Herbert D. Landahl (personal communication), is somewhat similar 
to the above but much more general. This model assumes that the 
A ‘triggers’? a chemical reaction which yields a great deal of 
energy. This may cause an increase of A’s velocity, but more 
important is the fact that the energy of reaction causes a local 
(with a radius of the order of angstrom units) heating of the region. 
This local heating may then cause a temporary ‘‘opening’’ of the 
membrane structure (i.e., by increasing the oscillation amplitude of 
the bonded particles). Thus, during this time A may pass through 
the opening. That is, the heating temporarily decreases the acti- 
vation energy necessary to cross the barrier. This model, as well 
as all others in this section, will be discussed at greater length 
in a later paper. 

Another model, which has been discussed under the general 
heading of carrier model, is a Maxwell Demon which causes a 
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transformation of A (such as an isomerization) that would permit - 
it to cross the potential barrier more readily. 

Two other models, similar to each other, have been suggested by 
Professor Joseph E. Mayer (personal communication). (A similar 
but much less detailed model has been presented by J. Danielli, 
loc. cit., mechanism e on p. 510). In these models the Maxwell 
Demon decreases the resistance of the membrane A by means of 
a chemical reaction. However, there is no increase of molecular- 
velocity and no local heating. In the first model suggested by J. 
Mayer, A reacts with the Demon which then allows A to cross the 
barrier. Later the Demon is restored to its original form by re- 
acting with a substrate. In a suggested example, the Demon is an 
integral part of the membrane and is located near the inner solu- 
tion. A forms a chelated complex with the Demon. This complex 
may then undergo a spontaneous rearrangement, in which a new 
bond is formed behind A, while a bond of the Demon molecule is 
broken in front of A. Then A is released. Since A cannot go back 
into the inner region, due to the formation of the new bond, it must 
therefore diffuse into the membrane toward the external region. 
When the Demon is restored to its initial state by reaction with 
some substrate, A is still prevented from going back to the inner 
region by the presence of the restored bond which was broken dur- 
ing the rearrangement and is now reformed. 


In J. Mayer’s second model A reacts with the Maxwell Demon to 
form a complex. Then, an enzyme reacts with this complex, al- 
lowing A to cross the barrier. Later a substrate molecule reacts 
with the Demon and enzyme and restores them to their original 
form. In a suggested example, A first forms a complex with the 
Demon (which is an integral part of the membrane), The enzyme 
then attaches itself to the Demon, in a manner such that A is sure 
rounded, and causes a rearrangement such that bonds are broken 
and a ‘‘hole’’ is opened in the Demon adjacent to A. The events 
following this are analogous to those of the model above. 


The Maxwell Demon is obviously very similar in action and in 
net result to the fixed carrier, although their detailed mechanisms 
are different. Both the Demon and the fixed carrier attach them- 
selves to A and then, without a translational movement ot the 
major portion of themselves, transfer A across a potential bar 
rier. In a given physical case the actual mechanism for active 
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transport will probably be a combination of a carrier and a Max- 
well Demon. FEowever, the Demon does possess the advantage 
(along with the fixed carrier) that its mechanism of affecting active 
transport need not involve the large diffusing molecules and/or the 
geometric complexities that most carrier models possess. 


The author wishes to express his gratitude to Drs. D. L. Cohn 
and R. Macey and particularly to Professor H. D. Landahl for his 
many helpful discussions and suggestions. 
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A learning curve derived by H. D. Landahl (Bull. Math. Biophysics, 
3, 71-77, 1941) from postulated neurological structures is shown to be 
derivable from simplified assumptions by introducing the information 
measure of the uncertainty of response. The possible significance of 
this approach to learning theory is discussed. 


H. D. Landahl (1941) derived a learning curve which relates the 
number of errors to the number of stimulus presentations in a mem- 
otization task, in which each of M stimuli is to be associated 
respectively with one of N responses: 


N 


Wa) aD va 


: l 
w= kb og 
where w is the number of errors, ¢ the number of stimulus presenta- 
tions, and kd is a single free parameter, presumably depending 
both on the characteristics of the subject and on M, the number 
of stimuli to be associated. 
By further assuming the relation 


k b = ne~ oM, (2) 


where now 7 and ¢ are two free parameters, I.andahl was able to 
fit three sets of data (for three different values of N =M) with 
only the two degrees of freedom. Each set contains 15 experi- 
mentally determined values (each based on 32-48 responses), to 
which the fit of the theoretical equation (1) with the right side of 
(2) substituted for 46 is practically perfect. 

The derivation of (1) is based on the properties of a postulated 
neural net. The postulates are, in effect, as follows: 


817 
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1. Only one of the set of responses is possible at each stimulus 
presentation.* 

2. There accumulates at the connection leading to the correct 
response an excitation factor, accumulation being linear with the 
number of correct responses. 

8. When the sum of the excitation factor and the net excitatory 
stimulation (due to the external stimulus at a connection) exceeds 
a certain (fluctuating) threshold, the appropriate response occurs. 

4. The threshold at each connection is a random variable with 
a Laplacian frequency distribution. 

5. There accumulates at each connection leading to a wrong 
response an inhibitory factor, accumulation being linear with the 
number of wrong responses of that type. 

Landahl averaged the inhibitory factors of all the N-1 ‘‘wrong” 
connections into a single inhibitory factor of ‘‘the’’ wrong re- 
sponse. Since, in comparing the theoretical curve with the data, 
the rate of accumulation of the inhibitory factor was set equal to 
zero, assumption 5 was not relevant for this case, and its as- 
sociated parameter does not appear in equation (1). This equation 
is based on assumptions 1-4, and the only relevant parameters 
are b, the rate of accumulation of the excitatory factor at the ‘‘cor- 
rect’? connection, and %, the parameter of the Laplacian distribu- 
tion; and even these appear only as the product bk, hence, as a 
single parameter. 

It will now be shown that equation (1) is derivable from the fol- 
lowing two assumptions. 

1, At each stimulus presentation, the range of responses funce 
tionally open to the subject varies, but the probability of each res- 
Sponse is the same and the ‘‘correct’? response is always included 
in that range. 

2. The amount of uncertainty of response (Shannon’s information) 
decreases proportionately to the rate of accumulation of erorrs. 

According to our assumption 2, we have 


-dH=kdwu. (3) 


*Strictly speaking, this assumption is a consequence of the experi- 
mental situation in which the subjects are directed to choose single 
responses. However, underlying this condition is a neurological model 
which makes such a situation inevitable. Therefore, we list this as a 
Separate postulate. 
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According to assumption 1, we have* 
H =-logP.,, (4) 


where P, is the probability of the wrong response. From (3), 
under the initial conditions H(0) = log N, we obtain 


H =logN -kw. (5) 
From (4), we obtain 
de 
Ne eA 8 ese (6) 


where c =¢ —w is the number of correct responses. Hence, taking 
exponentials 


d 
N — gk o ekt-ke. (7) 
Ne** de =e** dt, (8) 


Integrating (8), we have 
Nekc/k +A =e**/k, (9) 


where the constant of integration A is determined by the initial 
condition, c(0) = 0, to be (1 ~N)/k. Remembering that c =t-w, 
we arrive after re-arrangements from (9) to 


Now! ted (Ne= Ae Riri, (10) 
or, taking logarithms and solving for w, 


N 


shal. Ds 11 
(N-1)e"FF4 1 @ 


a 
w=-lo 
ke & 


which is exactly Landahl’s equation. Our parameter corresponds 
to Landahl’s single free parameter kd. If we make the same addi- 
tional assumption as Landahl, namely, & = ne~=™, we get Lan- 


*For convenience of calculation we measure information in ‘‘nits’’ 
rather than bits, which allows the use of natural logarithms. 


320 ANATOL RAPOPORT 


dahl’s two parameter equation which fits his experimental data so 
well. 

The introduction of the information measure seems to be parsi- 
monious in that the derivation appears to be independent of the 
several more or less ad hoc assumptions inherent in the neurologi- 
cal model, particularly of the rather specific Laplacian distribu- 
tion of threshold fluctuations. On the other hand, the present 
formulation may be only accidentally and formalistically related to 
the information measure H. It can be pointed out that we too have 
made an artificial assumption, namely, that the learning process 
consists merely of the elimination of wrong alternatives in an all- 
or-none fashion instead of (as in Landahl’s formulation) a gradual 
reinforcement of the tendency for the correct response at the exe 
pense of the tendencies for other responses. 

It can, however, be argued that in the case of memorization by 
human subjects this assumption is not implausible. It is equiva- 
lent to saying that once an association has been made, it remains 
for the duration of the process, so that (1) decrease of uncertainty 
comes only from errors, and (2) the only average effect of errors is 
to eliminate alternatives without biasing the remaining ones. 

If it should turn out that the information-theoretical model is a 
useful one in learning theory, the question naturally arises about 
the construction of information-theoretical models without the 
imposition of restraints on the relations among the probabilities of 
the responses. Generalizing equation (3), we would have, then, in 
the case of an information-theoretical model, an equation of the 
type: 

dH=f(pi, Payee, Pn-1, t) dt, (12) 
where the p; are supposedly independent probabilities of the re- 
sponses (p, is omitted, being determined by the rest). Since H is 


also a function of the p,, it can be omitted from the argument of f 
without loss of generality. However, d H is an explicit differential 


expression in the p,, namely, since Lid p; =0, 
dH =~ Dilog p;d pi. (13) 


Hence after substituting (13), equation (12) becomes a total dif- 
ferential equation in n independent variables. It has an integral 
only if certain stringent integrability conditions are satisfied. 
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When these are satisfied, the integral is an (n — 1)-dimensional 
manifold family 


ED GAD 95 sss Pais t= he (14) 


Such a manifold is a much ‘‘weaker’’ description of the learning 


process than is given, for example, by specifying all the proba- 
bilities as functions of ¢, 


Pi=p;(t) (i= 1,...,n-1). (15) 


Equation (15) determines a curve in n-dimensional space, that is, 
much more than an (n — 1)-dimensional hypersurface determined by 
(14). 

In our case we were able to derive a curve precisely because we 
put a strong restriction on the relations among the p;, according to 
which some vanished while the others remained equal. It may turn 
out, however, that in the case of very many possible responses 
(as, for example, combinations of muscular contractions in the 
acquisition of a muscular skill), no such specification can be 
made, but it may be possible to make a weaker specification, for 
example, to determine some function of the several responses 
(of which H is a special case) which does have a regular behavior. 
One feels that the rate of change of H is related to the rate of 
‘¢information storage’? in the organism. If no more than a many- 
dimensional manifold can be determined as characteristic of a 
given type of learning process, it is futile to look for learning 
curves. All that can be determined is a learning ‘‘hypersurface,” 
any one curve along which may be the learning curve pertaining to 
an individual or to a particular realization of the learning situation. 
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In connection with previous studies (Rashevsky, Mathematical Biology 
of Social Behavior, chap. xii), a situation is investigated in which the 
two mutually exclusive possible behaviors of a society consist of the 
desire to keep the present socioeconomical situation and the desire to 
change it in any way. The psychophysiological tendency ¢ towards 
either of the bebaviors is considered to be proportional to the difference 
between the actual income 7 of the individual and his needs 7% As- 
suming that the distribution function N,(@’) of the needs is a given 
characteristic of the population, it is shown that the distribution func- 
tion M(¢) of ¢ in the society can be derived from N,(¢’ ) and from the 
distribution N,(z) of the incomes 7. A particular case is worked out as 
an example. Conditions of stability of a socioeconomic structure are 
studied in their dependence on the income distribution, 


In previous publications (Rashevsky, 1949; 1951, chap. xii, re- 
ferred to as loc. cit.) we have discussed the effects of imitation on 
behavior of large groups of individuals. We considered the situa- 
tion in which each individual is characterized by a tendency ¢ 
towards one of two mutually exclusive behaviors. This tendency 
is measured by the difference in the excitation of two centers 
which mutually inhibit each other. The excitation of one center 
produces one of the two mutually exclusive behaviors; the excita- 
tion of the other center produces the opposite behavior. When both 
excitations are equal and ¢ =0, the individual has no preference 
towards either of the behaviors. If ¢>0, the preference is for one 
behavior which we may denote behavior [; if ¢ <0, the preference 
is for the other behavior which we shall denote as behavior Il. For 
‘a given ¢, an individual exhibits both behaviors with a relative 
which is a function of ¢, according to the well-known 


c 
oe. 1938; Householder and Landahl, 


Landahl’s equation (Landahl, 
323 
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1945). If f =0, the relative frequencies are equal: the probability 
of each behavior is 4%. As ¢ increases, the relative frequency 
of behavior I increases, that of behavior II decreases. As ¢ de- 
creases, the reverse takes place. 

In loc. cit. we considered a population in which the quantity ¢ 
is distributed at random according to some distribution function. 
We assumed that the distribution function is symmetric with re- 
spect to ¢ =0, having a maximum for ¢ =0. We considered the 
distribution function given by the following expressions: 

If N(¢) d(¢) denotes the fraction of individuals with a d between 
¢ and d +d dq, then 


For ¢ > 0: Ng) = Se 74; (1) 


For 6 <0: N(d¢) = ze od, (2) 


Subsequently H. G. Landau (1950) considered the case of a 
normal distribution and even generalized our results to any distri- 
bution functions. 

The purpose of this paper is to study a distribution function 
which is largely determined by economic conditions. 

Let an individual have an income i. If z is not too small, then 
some individuals may find it greater than their minimum need, some 
may find it just equal to their minimum need, and others may find 
it less than what they need. Let us denote by i’ the minimal need 
of an individual, in terms of monetary income. There will be in 
general very few individuals whose ¢’is either unusually small or 
unusually large. The majority of individuals will have some aver- 
age 7’, Hence, we may expect that in general the distribution func- 
tion N,(2’) will be equal to zero for 7’.= 0; will have a maximum for 
some value 7’,.and then drop asymptotically to zero with increasing 
us 

Such a situation is, of course, highly oversimplified. We assume 
that the minimal need of an individual is an endogeneous character- 
istic of the individual. Actually it is more than likely that the 
minimal need, as it is felt subjectively, will depend on the incomes 
of other individuals with whom the given individual comes into 
contact. An individual who sees another with a much greater in- 
come than his own may consider his own income inadequate. The 
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same individual living among poorer individuals may consider his 
income as adequate. 

Such considerations lead to problems of considerably greater 
complexity than the one discussed here. We will therefore limit 
ourselves here to the simpler case in which N,(2’) is a character- 
istic of the individual only. 

It is not implausible psychologically to assume that an individ- 
ual, for whom i =2%4, will be satisfied with the existing economic 
state of affairs which leads to this equality of i and 17 However, 
while not dissatisfied, he may have no particular reasons to strongly 
like the existing economic state of affairs. An individual with 
t <2’ will have definite reasons for being dissatisfied with exist- 
ing conditions, while an individual with <> 7’ will feel the more 
strongly in favor of the existing conditions the stronger the in- 
equality 77>. 

The efforts to maintain the economic status quo and the efforts 
to change it in any way may be considered as two mutually ex- 
clusive behaviors which we will denote correspondingly by I and 
II. The above considerations then suggest that as the measure of 
the reference towards behavior I or behavior II we may conven- 
iently choose the quantity 


gb =i-i’ (3) 


We will now show that the distribution function N(¢) can be 
derived from the distribution function N,(i’) and the distribution 
function N,(2) of incomes. 

Let us first consider the somewhat more general case in which 
¢ and 7’ are correlated. Denote by G(i,z’)dide’ the fraction of in- 
dividuals who have an income between 7 and ¢ + dz and a need be- 
tween 2’ and 2’+ dz’ 

We then have 


[ ecinai = NG, (4) 
0 
and 


| G(i,i’)di’ = No(i). (5) 
0 


326 N. RASHEVSKY 


Because of (3) the number of individuals who have an 7 between 
¢ and z + di and an 2 between 2’ and 2’ + di’ may also be said to be 
the number of individuals who have an 7 between i and i +di and 
a ¢@ between ¢ and ¢ +d¢d, where ¢ =i —7% 

Now we must distinguish between the cases 4 >0, and ¢ <0. 

When 7’<z and therefore ¢ >0, we find for the fraction G,(z,¢) 
didd of individuals with 2 between i +di and ¢ between $6+dd¢ 
the expression 


G ,(i,@)didd = G(i,i — d)didd. (6) 


To obtain the distribution function N(d) we must integrate 
G ,(i,¢) with respect to z, from i = ¢ toi = «. The lower limit must 
be ¢=¢ and not i = 0 because, for a given d >0 and a given Zz, 
the quantity 7 ~ ¢ in (6) cannot be negative since it represents the 
need 7% 

However, the same result is obtained obviously by considering 
the fraction Go(¢,2’) of individuals with a given 4 and given 7’. We 
find then 


Go( bt )dddi’ = G(i’ + 6, dddi* (7) 
Now we have, 
for ¢ >0: N(d) = fows ot )di % (8) 
0 


When z’>z and ¢ <0 then it is more convenient to use expres- 
sion (6) and we find 


eo 


ford <0: M(¢) = | Gi aha (9) 


0 


For the special case in which i and 7’ are not correlated and 
which we shall consider here in more detail, we have 


G(i,2’) = 1(¢)N o(7). (10) 


According to (8) and (9) we have: 
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For $>0: NS) = Wald) = fv s(iQnalirs gai’ (11) 
0 


co 


For ¢<0: — N(¢) = N_(¢) = i N (i - AN, (di. (12) 


0 


If, as we said above, we consider the need 7’ as an inherent 
property of the individual, then N;(2’) will not change with chang- 
ing economic conditions, but N,(z) will. Therefore, as N,.(2) 
changes, ‘the shape of the function N(¢) will also change. But as 
we have seen in loc. cit. some changes in N(d) may result in 
sudden changes of the socioeconomic behavior. 

To investigate those situations under specific conditions, we 
shall choose some special functions for N1(2’) and N.(z). We shall 


arbitrarily take 
oy , 
N,1(’) oe 2 en bi e (13) 


It satisfies the condition 


[ N1(i)di’ = 1. (14) 
0 


As we have seen elsewhere (Rashevsky, 1951, chap. ix) the 
income distribution may roughly be represented by a simple ex- 
ponential. A better representation is given (Rashevsky, 1950) by 


F(a) = Cle“ XO ~ [14 B(i- Ale PMI, (15) 


where C, o, B, and A are constants. In particular A denotes the 


smallest subsistence amount. 
For simplicity we shall, however, use a simple exponential and 


put 
N,(i) = ce ~°', (16) 


which also satisfies the condition 


60 


| N2(i)di = 1. (17) 
0 
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Introducing (13) and (16) into (11) and (12) we find (Peirce, 
1929, page 63, formula 493): 


8 é 
For ¢ > 0: N(#) = N,.(¢) “ire eb (18) 
b8 1 ) ¢? 
For d mes N(¢) = N_(¢) 7 (eno haat + = eo, (19) 


It can be shown that, as of course should be the case, the area 
under the course N(¢) is equal to 1. In fact the area under the 
branch N_(¢) is equal to 


i : bc be c (20) 
= ——_ + ——__ 4 —__, 
a AOS aay eRe ieee 
The area under the branch NV, (¢) is equal to 
; wt (21) 
N dd = ——___., 
The sum of (20) and (21) is equal to 
ee b2¢ be () 
a pe lk 
(b +c)’ ae +ce)8 ar +c)? bre 
(22) 
6° 4570 +bce(b +c) +e(b+c)? 52 4 352c +302 +3 
(b +c)§ (b +c)8 oe 


Though the distribution function N() is composed of two dif- 
ferent analytical functions, it is continuous everywhere and has 
a continuous first derivative. In fact, we have from (18) and (19): 


63 
N(0) = N (0) = N_(0) = Rae (23) 
We also have 
dN. (4) beet 
damascene (24) 
aN_(d) 8c ~c cd om rv 
meee 
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From (24) and (25) we find 


EO) {2x0} {Ee 5 end: Oo ag 
ES be dd el de Se (26) 


The function N4(¢) is monotonically decreasing between ¢ = 0 
and ¢ =o, The function N-(¢) has a maximum. As a matter of 
fact N_(¢), when extended into the region of positive values of ¢, 
also has a minimum. This can he seen by equating the right side 
of (25) to zero. 

This leads to 


c 
jp = —————— = 0. 27 
$ (b +c)? (27) 

The first and last coefficients of the quadratic equation (27) are 
of different signs. As can be readily verified the positive root cor- 
responds to a minimum, the negative to a maximum. For c =0 
there is only one root, namely, 6 =0. For c = « equation (27) 


becomes 
5b? +6 =0, (28) 
the roots of which are 
¢, =9 and ¢, = =. (29) 


The full line in Figure 1 represents the general shape of the 
function N(¢). The broken line represents the extension of N_(¢) 


FIGURE i. 
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into the region of positive values of ¢. The curve is drawn for 
the casec <b. 

The parameter 5 reflects the distribution of the individual needs 
and, with the simplifying assumption made above, should be con- 
sidered as a constant for a given population. On the other hand, 
the parameter c which reflects the distribution of incomes and 
whose inverse value represents the average per capita income 7, 


=ary 
t= 


Cc 


(30) 


’ 


may vary with changing socioeconomic conditions. We shall there- 
for study the effects on M(¢) of varying c when d is kept constant. 

From the discussion above it follows that the function M(d¢), 
being composed of two analytically different branches, is never 
symmetric with respect to any value of ¢, regardless of the choice 
of the parameters 6 and c. It is, however, possible to choose such 
a value c* of c that the areas under the curves Ni(¢) and N_(¢) 
are equal, To that end we equate the right sides of (20) and (21). 
We thus obtain: 


b3 b?c be c 
(CPanel ae be 
or, clearing fractions and collecting terms: 
c® +8 7c +3 bc? =53, (32) 
The equation (32) may be written: 
(6+c)? ~53 =53, 
which gives 
ada 
(b Fe)? 9" 
This gives for the value c*: 
o* = (¥/2 = 1)b = 0.26. (33) 


We now have the following general picture of the behavior of 
N(¢) with varying c. 

For c = c* = 0,265 the mean value of @ is zero. As c decreases, 
tending to zero, the area under the curve for positive values of ¢ 
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increases as compared to that for negative values. From equation 
(31) we see that c can always be chosen sufficiently small so as 
to make the ratio of the left side to that of the right as large as we 
wish. This means that c can always be chosen small enough to 
make the area subtended by N,(¢) as close to 1, the total area 
subtended by N(¢), as we wish. Because the total area remains 
equal to 1, the ordinates of N(¢) tend to zero as c decreases. At 
the same time, as seen from the discussion immediately following 
equation (27), the position of the maximum of N(¢) tends to ¢ = 0. 
The situation for rather small values of c is illustrated in Figure 
2. 

As c increases, the branch N,(¢) shrinks to zero, and N(¢) tends 
to N_(¢). The function N(¢) has now, as seen from (29), a mini- 
mum at ¢ =0 anda maximum at 6 =—2/b. All of the area subtended 
by N(d) tends now to be under N-(¢). The situation is illustrated 
in Figure 3. The three figures are not made exactly to scale. 

We now may use the function N(¢) for the discussion of imitative 
behavior in the same manner as we did previously (Rashevsky, 


1951, chap. xii). 
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The increment y of d, due to imitation, is again given by equa- 
tion (8) of loc. cit. which we shall rewrite here again. If X and Y 
denote the numbers of individuals exhibiting behavior I and I], and 
if A and a denote neurophysiological constants, we have: 


OY ieee 34 
A(X -Y) -ay. (34) 


Denoting by x(¢)d¢ the number of individuals with a given ¢ 
which exhibit behavior I and by y(d)d¢ the number of individuals 
which exhibit behavior II, we may again, like in loc. cit., write 


+00 
X-Y =N, [xo [2P (4,9) - 1 d¢, (35) 
where 
P, (¢,¥) =1~Y%e “*F + for g>-y, (36) 
and 
P (ps) = he *F * for 6 <-y, (37) 


and NV, is the total number of individuals.* 

The integral in equation (35) can be readily evaluated in a simi- 
lar manner as has been done in loc. cit. However, the final expres- 
sions are somewhat unwieldy and have little interest at present. 
We shall therefore study the solution of (34) without actually 
evaluating the integral in (35), for the particular case when & in 


(36) and (37) is very large. In that case we have with good ap- 
proximation 


P(d,y)=1 for $+y>0; (38) 
P.(¢,¥)=0 for d+y<0. (39) 


The meaning of (38) and (39) is that every individual with 
@+Y>0 exhibits behavior I; every individual with @é+w<0 
exhibits behavior II. 


*In loc. cit. N(¢)d¢ referred to the number of individuals with a given 
$; here it refers to the fraction of individuals. Therefore in equation 
(20) of loc. cit. from which we obtain our equation (35), N_ is absorbed 
in (4) [cf. equations (11) and (12) of loc. cit.]. 
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With the assumption (38)-(39) about the P, function let us return 

to the distribution N(¢) as given by (18) and (19) in the absence 

of imitation. The difference X —Y is then simply equal to the dif- 
ference of the areas 


oo 


i N($) a4, (40) 
0 
and 
O 
[ nw ae, (41) 
multiplied by N,. Imitation has the effect of shifting the whole 
N(q¢) curve by an amount 7% parallel to the ¢-axis. The value of 
A(X —Y) for any value of w is now given by 
oo — 
A(X ~ Y) =F(¥) - aed cen - [os (42) 
~y wes 

From (42) it follows that the derivative dF/dy is equal at each 
point 2AN ,N(¢). 

Let us now consider the case represented on Figure 2, when 
ce <<b. Then F(0) is almost equal to +1. As yw increases from 
zero, F(w~) approaches 1 asymptotically. We have F(y) =0 for 
some rather large negative value of y. For more negative values 


FIGure 4. 
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of w, F(w) is negative and tends asymptotically to —1, but the 
asymptotic approach to ~1 is much slower than the asymptotic 
approach of the positive branch of F(w) to +1, because N(d) de- 
creases much more rapidly for negative than for positive values of 
Ww. The F(W) curve in this case looks like the curve shown in 
Figure 4, The highest value of dF/dy = 2AN,N(w) will be for a 
slightly positive ~, where F (7) is already almost equal to 1. Thus 
the inflection point is very near to the asymptote +1. 

The equilibrium values of yw are determined by the points of 
intersection of F() with the straight line ay. Except for une 
usually small values of a, only one point of intersection, 01, is 
present, It corresponds to a stable equilibrium in which the vast 
majority exhibits behavior I; in other words, this majority desires 
the preservation of the status quo. For exceedingly small values 
of a, three equilibria are possible of which one is unstable and 
two stable. Such is the situation for a small c, that is, according 
to (30), for a large enough average per capita income 7. 

Let now c increase and tend to infinity. From (18) and (19) we 
see that the function N(¢) tends to 


p3 
ee 6%: 4 <0; (48) 


which is represented in Figure 8. Now we have in the limit F(0) = 
—1. The function F(¥) is equal to zero for a positive value w* of 
¥. The maximum of the derivative dF/dy, that is, the inflection 
point, corresponds to a positive value of ¥, such that 0 < w<*. 
We have a situation which is represented in Figure 5. 

If we have originally a small c so that only one equilibrium 
point O,, corresponding to behavior I, exists (Figure 4) and if c 
gradually increases, then eventually the situation represented in 
Figure 5 will be reached, in which, for not too small values of a, 
only the equilibrium point O, exists which corresponds to the 
majority exhibiting behavior II, At some intermediate stages a 
situation will exist in which the positive branch of F(y) is just 
tangent to the straight line ay and at which a sudden change from 
behavior I to behavior II will occur. 

It must be noted that, whereas the root of F(w) =0 tends to 
~ee as c tends to zero, that root tends to a finite positive value as 
¢ tends to infinity. Therefore if a is sufficiently small, the straight 
line ay will never cease to intersect the positive branch of y, as 
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FIGURE 5. 


shown by the broken line of Figure 5. We now have three points of 
equilibrium, 0;, O02, and O3, not shown in the figure, of which 01 
and O, are stable, O, unstable. In that case, if we originally had 
the situation represented in Figure 4, then no matter how large c 
becomes, the point O, will always remain stable. A large c, as we 
have seen, means a small per capita income 7% and a concentration 
of wealth in a small fraction of the population. We thus come to 
the following result. If a is sufficiently large, then for large per 
capita income the population will be satisfied with the existing 
socioeconomic status. But as the per capita income drops below a 
threshold value, a general discontent and opposition to the exist- 
ing order will suddenly occur. If, however, a is sufficiently small, 
then no matter how low the per capita income drops, the majority 
will remain satisfied with the existing state of affairs. A ‘‘revolu- 
tion’? will never occur for sufficiently small values of a. 

Because of the oversimplifications in our basic assumptions, 
those conclusions cannot be directly applied to any actually known. 
situations. Their interest, however, lies in the circumstance that 
we have here the effect of a psychophysical factor, a, upon the 
stability of an unsatisfactory situation. 


This investigation was aided by a grant from the Dr. Walter OF 
and Clara A. Abbott Memorial Fund of The University of Chicago. 
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A method is developed, which is essentially an elaboration of the 
Pitts-McCulloch net, to give a systematic and effective account of be- 
havior in mathematical terms. In particular, in this paper an interpreta- 
tion is given cf an expectancy-type theory of behavior with suggestions 
for its extension. 


This paper aims to give a method for developing mathematical 
logical networks and showing the manner in which they can be used 
in interpreting a particular central fragment of a theory of behavior. 

Notation. An excitation input fibre ending on an element has a 
filled-in triangle drawn inside the large circle which represents the 
element, while an inhibition input fibre ending on an element has a 
small open circle drawn inside the larger circle. The threshold 
number, which states the number of excitatory inputs that have to 
fire simultaneously to fire the element, is given by some real inte- 
eral number n written inside the large element-circle. The output 
fibre has no special ending inside the element (Stewart, 1956). All 
the churacteristics named so far are incorporated in Figure 1. 


FIGURE 1. 


Every such network that can be drawn has an accompanying 
mathematical logical formula which the network is said to realize. 
Thus the formula that gives the condition for F to fire in the net- 


work of Figure 1 is: 


2 ta 
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Certain conventions will now be briefly listed that apply to all 
such networks: (1) We shall assume that all elements of a network 
are either in a state of excitation or they are not. (2) Errors will 
be assumed not to occur. It being noticed that J. von Neumann 
(1952) has made it clear how errors could be introduced into the 
system or, more accurately, the methods by which we can ensure 
efficient working with inefficient components. (3) Time will be as- 
sumed to be in equal units which will be taken as the firing time of 
any of the elements. (4) We shall talk of inner and input elements 
according to whether or not they have input wires. (5) The same 
output wire may be connected to several inputs, but any one input 
is assumed to be connected to at most one output. 

Logical Constraints. The whole of the present discussion takes 
place within the framework of S. Kleene’s (1951) regular nets. We 
shall follow Kleene’s definition of a regular event and say that an 
event is regular if and only if it is made up of definite events E,F,... 
by use of the operations...v---, ...--- and ...*---, where 
in terms of — and F the operations are written E v F, EF, and 
E*F, ‘‘y’’ is the Boolean ‘‘or,’’ and EF means that E occurs having 
been immediately preceded by F. E*F is EE...EF. 

All such regular events are primitive recursive and realizable, 
and a definite event is any event that is representable by the firing 
or nonefiring of an inner neuron at some time in the specified past, 
or present, of the organism. Thus serious logical difficulties are 
avoided by keeping within the domain of regular events. 

The representation used will thus be in terms of the proposi- 
tional calculus, with a fragment of the lower functional calculus, 
all subscripted for time. 


Before proceeding to generate theorems, or organs, on the above 
basis, let us outline the two postulates of behavior that are inde- 
pendent developments of psychological theory and that have been 
written up in detail elsewhere (George and Handlon, 1955, 1956). 
Let us start with some glossary definitions (George, 1953). 


Behavior Postulates. D1. A belief (B)isa relatively permanent 
state of the organism which represents the association, through ex- 


perience, of the activating stimulus state (S,) with a means-oute 
come (Rk, —S,). 


D2. An activating stimulus state (S,) is compounded of all the 
effective stimuli, motivationally neutral and otherwise, that ef- 
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fectively change the state of the organism. These substimuli may 
be internal or external to the organism. 

D3. Motivation (M) is a complex state of the organism that is 
necessary to a change of behavior. It manifests both primary and 
secondary drives and operates as part of S,. 

D4. A stimulus complex (S*) is that part of the S, that is not 
contributed by M. 

D5. An expectancy (E(R, —S,)) or more briefly (E), is a rel- 
atively temporary state of the organism, elicited by S,, that has 
the same content as its correlative B. It is a B-in-action. 

We need not bother too much here about a means-outcome except 
to see it as a state of the organism that is contingent upon the 
action of the organism. A means-outcome is a combination of an 
action and the result of that action. The process of sensitization 
is one that refers to the beliefs that are immediately elicited by 
any S;. 

Now let us state the two main postulates: 

P1. A necessary and sufficient condition for the occurrence of 
an expectancy in an organism is that for M > & and S* > 1, where & 
and 1 are constants dependent on the circumstances, there is at 
least one B, activated by the S;. The Sj, in general, translates a 
range of B’s into a range of £’s, Thus after sensitization there 
exists a partially ordered set of E’s. 


P2. The particular range of B’s translated into H’s are acti- 
vated as a ratio of recently successful responses n, to total num- 
ber of past responses ny. Thus in any interval (¢),¢,), the strength 
of the belief (SB) can be written 


SB; (¢1) = SB; (to) +) SB: (n,t1), 


where the range of z is over all the sensitized B’s in time interval 
(to,#,) and n =n,*/n,*, where n,* and ny* are values of n; and ng 
in the same interval. It should be clear that any stimulus that 
arouses a set of B’s will lead to a response, or set of responses 
if the B’s aroused are consistent with such responses. Reflective 
behavior involved further stimulation where there is an awareness 
of inconsistencies in the B’s aroused. 

Logical Networks. Our next task is to develop further the organs 
of the logical networks. Let us describe what we shall call a B- 
network. Figure 2 shows the simplest kind of such general net- 
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work, where the relevant formulae are 


ia Si sae (2) 
L} = Se ae Shae : (3) 
[ios ene (4) 
Mie~ Lie ~L? + (Ki, vat) Me); (5) 
Rises en: (6) 
Re tan Stee, (7) 


R1? a R21 a (Se : M2) vy (See 5 Se Vv Si Mi 3). (8) 


FIGURE 2, 


By substitution of (2), (3), and (4) in (5), and by further substitu- 
tion of (5), (6), and (7) in (8) we get: 


ibis = (Si bear Loe athe Lee ; 
vise v (47*) M} (t*)) v (Se) : S2,) Vv Sap .~ Des . (9) 
~ Lt 4) +(Ki_g v (At*) M? (t*))), 


where ¢, < ¢* < ¢, and ¢, is the initial time for the machine’s opera- 
tion, and there will be no firing of L’ and L? between ¢* and ¢. Or 
t*represents the last time K' fired and this is retained in the loop 
of M}, 

The formulae (2)-(9) define a B-net which has an associated 
B-function B(n,m) where n denotes the number of simultaneous 
firings of S' and S? that are capable of being recorded, and m 
denotes the number of firings of either S! and S2 alone that are 
capable of being recorded. Thus B(n,m) describes the limits of the 
double-counter of conjunctions and disjunctions of the two stimuli 
or inputs, 
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The simplest sort of B-function associated with the Benet of 
formulae (2)-(9) is of the form B(1,0). (See Fig. 2). Figure 3 is a 
realization of B(3,3). 

It is easy to see how the generalizations of these counters 
should take place in terms of logical networks. This implies the 
addition of one formula for each element that is added to the nete 
work, 

The number of elements required for the realization of B-nets 
capable of recording greater numbers of conjunctions and disjunc- 
tions is easily achieved. Thus: 


B(1, 0) =9 
B(2, 0) = 10 
B(3,0) = 11 
B(n,0) = n+ 8 
B(1, 1) = 10 
B(1, 2) =11 


B(1,m) =m +9 


B(n,m) =n +m+8 


So we have a simple general formula for the number of elements 
needed on this method of construction. Thus it is easy to compute 
the number of elements needed to realize a set of B-functions, 
provided we know the values for n and m of each function. This 
so far only holds for B-functions connecting 2 inputs S$ and Si, 
(2 4 7); we shall consider more general B-functions shortly. 

The effectiveness of a Benet is given simply by a single number 
r where 

= (n ~ m). 

Thus B,(r) is the single-valued function, having as its possible 
domain all real integers, that describes the state of any B-unit at 
any time ¢, 

A special modification of the Benet should now be mentioned 
where the original conditions are changed. This can be shown by 
substituting for (8), a new formula: 


Byte Si -1 ¥ Sto, Mig) ¥ (SE pastor y v:(SieMioy, ome 


fori = 1, j = 2. 
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This means that stimulation of S! will always fire R*/ as well 
as f/, This particular B-net is therefore the analogue of the con- 
ditional response or any relation, such as that a new stimulus has 
come “‘to stand for’? an existent permanent S-R connection. 

We can, using the above B-nets, now develop a conditional prob- 
ability machine of the kind suggested by A. Uttley (1954, 1955). 
Thus consider any general B-function (9) with values n =2 and 
m =y, and let the occurrences of S! and S? be given by the follow- 
ing table: 

Oot li. twee Lae cad ek: 
eet eee 1s ei 

In the above table, in the 11th column, for example, if oe fired, 
the probability of s’ firing will be given by some number p, which 
may be calculated according to the particular nature of the B-net 
and for particular choice of z and y. But in particular for the most 
general form of Benet, p will be 6/10 and S_ will fire if 2 ¢ 2 and 
y need take no value other than zero. If the above table had 
started: 

Seo USO Ue ta tte) 

Seles Wont hele Re aR coi 
then S!’s firing at the 7th column will not cause S? to fire, since 
the non-contiguous firing of S' and S? up to the fifth occurrence 
did influence the state of the network, provided y took some value. 
y = 4 would be sufficient to keep the full record of non-contiguous 
firing. 

Classification and Control Machines. Let us define a classifica- 
tion and control machine of n-inputs as that collection of B-nets 
which connect each of the n-inputs in every possible combination 
1,2,..., 7 at a time. 

Thus a 2-input machine will be a simple B-net given by (9); a 3- 
input machine will have 


gC 1 +3Ca + 9Cg = ¢ direct outputs, 


where three of these come straight from the inputs, (i.e., Cal he 
R®), and the other four from the four B-nets connected to the three 
9 


inputs S!, S?, S? in pairs and in a triad. 
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We can represent a 3-input machine by the diagram in Figure 4, 
where the squares are taken to represent the B-functions connect- 
ing inputs as stated in the relevant squares. 


R25 
R 


FIGURE 4. 


The words ‘‘direct output’? refer to the immediate results of 
classification on inputs. Clearly the direct outputs can be mapped 
onto any number of actual outputs, although this number will gen- 
erally not be less than the number of inputs. For our purpose we 
need not here consider this problem and may therefore ignore any 
motor classification system. 

In general an n-input classification and control machine will 
have 


nC 1 t nU'2 teeoot now 
direct outputs and thus will have 
nC'4 =f uaa t+eeot aun —-n 


Benets, each represented by a closed square. 

Although the Benets so far discussed relate to 2 inputs, the 
generalization to 3,4,...,S inputs is quite simple. The only 
elements affected are the L-elements. One extra L-element being 
needed for each extra input. In general there are two such L-ele- 
ments L,, L, in all B-unit networks combining 2 inputs. They are 
the recorders for the firings of inputs independently (see Figs. 2 
and 3), 

One important thing to be noticed is that motivation is not cas 
tered for by this ‘‘paper machine” (Turing, 1955) as yet. We can 
of course easily introduce motivators and such would be necessary 
to bring the model into line with the facts of behavior. Actual 


motivators are dependent on the materials of construction of actual 
machines, 
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The above point is important because it accounts for the fact 
that all contiguous events are not retained by the organism. The 
ones retained are those that are reinforced in terms of the motiva- 
tions of the organism. Thus motivational stimuli can, for example, 
be thought of as S_’s firing into all the B-units and into the K- 
elements of the general B-unit and changing the threshold of these 
elements by one apiece. 

The process of reinforcement will operate by firing more M-units 
with each contiguous occurrence of S’s together and diminishing by 
some multiple all S’s that are not fired contiguously with certain 
S’s that are fired. 

This operation of inhibiting connections not reinforced is per- 
formed automatically when B-units are connected as they are in our 
classifying-and-control machine. The point is that the mere occur- 
rence of two stimuli simultaneously would lead to the setting up of 
associations as mirrored by B-nets. This would be unrealistic 
from the behavioral point of view and thus only those associations 
that stop the motivator input lead to strong associations. Indeed 
at this stage this means that, in the machine, the response must 
inhibit the motivator and thus stop the input, bringing about a state 
of homeostasis in the machine. 

The formulae defining a B-function B(n,m) is given by: 


Ke Sera Sea! (11) 
Pease Si} (12) 
Lie So Sa (13) 
Rp = Spa43 (14) 
Re = S143 uy 


GPa (K LewGi) C1). OL ae hi Ok 08) 
0? = (KL, v(3t*) C10") + (34?) 707) (Ly Lis OD 


Om =((KL, Cnt) v (ae") CM) (Lt: mots (18) 


— 1 
n n 1 
Doereth tay G2 gh v (AP) B™” 
QP le yt CRS); 


(19) 
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D™ = ((L2., vb2.,) v(at7 st B) DP EEG? TR) (oak ease) 
Re = (Same *Cm1) Ys (Srey S254) Ms (Sty Cin) (22) 


It should be noted that there is the possibility of developing this 
theory in terms of matrix algebra, in a manner that has been pre- 
viously suggested in this journal (Landahl & Runge, 1946), Thus 
the complete ‘‘structural’? matrix F for B(n,m) is: 


Lite h +0... 050s On ce 0) Ont ene alae 
to 15 120-000... 022020 20. MLO ae 
0 One 08 Aha ee el tlre 0 0 
Oe 0 0) aca eS ee ke ah be Basal | ww 
0. 0 (: 0. —1 —laaIi w= b> defi ie. corre 
0-050 ot- ahie.0.5.0 5 20 20, Ors. .Oni Ou 
OS Os OF Oita lar. On Oo OO ar On IROMEO 
O20 27:00 Se Ona ra Oe Ola OM Ore nr On |e OmO 
OG RO Os Om Ose Oc toll Oe Obs. Ont Ome 
Ose OO Os OO Oe eee Orr Ons Om Oma 
OF. OO 0 0 Or. Ons Oe eer areas | POMEO mman 
OF 0% 05 08 20's -0- 7.00 OM lar Or tO 
Oo 00-1 


0 


In practice we shall be concerned with the matrix of the machine 
rather than each Benet. Thus in the simple 3-input machine, suc- 
cessive “‘status’? matrices will be summed to give the probability 
p of a particular output firing. Hence the first four status-matrices 
for 3 conjunctions S, +S, and 1 conjunction S,-S.-S, may read: 


The probability of firing p, in the present design, is such that if 
r>0, where r=(nm), the output fires. This means that p = 
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a/y +a, where n=a, and m= y, and that p must be > % for the 
output to fire. 

It is now fairly easy to see how the B-nets so far developed can 
be interpreted as the beliefs of P/ and P2. We set up the inter- 
pretation by saying that the M and the S* are the input fibres of 
input elements that are sufficient to fire the network thus eliciting 
a B, or range of B’s, and thus a range of E’s. The actual ordering 
of the B’s can be made on the basis of the conditional probability 
ratio of number of successful occurrences to total number of occur- 
rences. This, in general circuit terms, merely reauires that we 
have counters. In terms of the networks previously outlined it 
means that, if we suppose that initially all the @’s and y’s are 
zero, then the non-occurrence of, say, S; with S,; implies # is 
diminished, by one; whereas the contiguous occurrence of S$, and 
S, adds one to the y-number. Or, more simply, we may regard the 
y’s as negative values of x Thus there is but one number that 
now ranges over positive and negative numbers, and its value r = 
k, say, will be the weight of the contiguous firing over or under 
the total firing of a particular S; with respect to some S;. 

The above gives an exact measure of conditional probabilities 
based on experience. It, in fact, represents the factor usually re- 
ferred to in psychology as frequency. The factor of recency can 
be given weight by causing the effect of the early x-values to be 
diminished in time. If this were not so, then it would be virtually 
impossible to break anything that might be reasonably described 
as a habit. In practice this means that a temporal delay in the use 
of certain circuits must gradually and increasingly lead to stimuli 
of the x-reducing character to be fired rather as an alarm clock 
fires. This is, however, a matter of minor detail. 

After the existence of a partial set of E’s, the process is more 
complex, or net, according to whether the set is a unit set or a 
consistent set. If in fact only one sensitized B occurs, or a set 
of consistent B’s occurs, or to put it the other way around, if no 
inconsistent B occurs, then no further operation (‘‘choice’’) is nec- 


essary. This sort of situation possibly refers to most acts of the 


habit-type. : 
If, however, there exists a proper set E,, not all of which can be 


acted upon, then the process is more complicated and goes beyond 
the postulates. A process of evaluating occurs, which involves 
what we call ‘‘consciousness”’ in organisms, and probably involves 
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a further scanning and weighting process of our originally sensi- 
tized set of E’s. The fact that certain circumstances give rise to 
doubt about the means-end to be acted upon will be due to the fact 
that certain sets of E’s may be reinterpreted, or changes may occur 
in the classification system which is responsible for n,, n, and n, 
or there is perhaps no (classified) similar previous experience, or 
perhaps a lot of experience of roughly equal strength. This diffi- 
culty arises with a very complex organism and occasions what we 
call hesitation or doubt and is made very much more complex by 
the intrusion of consciousness which at the network level implies 
a further firing of S’s that may change the probability situation 
quite quickly, even as the counting is occurring. 

We cannot give any detailed picture in this note of this more 
complex process which mirrors reflective action, except to say 
that it involves a repeated scanning of the partial set of E’s in 
the same essential manner as the ordinary habit-eliciting process. 
Thus the only difference is one that should be capable of being 
catered for by viewing many B-nets connected in a large machine. 
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ERRATA 


1. To the paper by Eugene Ackerman entitled ‘‘Resonances of 
of Biological Cells at Audible Frequencies’’ (Vol. 13, pp. 93-106.) 
One term was omitted from equation (5) on page 96. The cor- 
rected equation should read: 
rk F) (. a Z oR 27 
in 9 36) * aF sin? 6 dy? a? 


~ a? sin 6 0 00 


Pye koee & 
This last term was also omitted from the last three equations on 


page 97 and the first four on page 98. 
The second equation on page 98 should read: 


‘ n Tn 
wo ( + 2 =— (n -1)(n + 2). 
fa n+1p,] p,@ 


The fourth equation on page 98 should read: 
w? = 4.8 Ta *. 


Although this increases the numerical values of T by a factor of 
three halves, it in no way alters the interpretation which was based 


on orders of magnitude only. 

9. To the paper by Eugene Ackerman entitled ‘*An Extension of 
the Theory of Resonances of Biological Cells: I. Effects of Vis 
cosity and Compressibility’’ (Vol. 16, pp. 141-50). 

This same error reappears again equation (14) on page 144, A 
term should be subtracted from the right hand side having the value 


. 2 ris 


This error was masked by a numerical error on page 148, where the 
simplified equation of motion is shown correctly as: 
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38. To the paper by Eugene Ackerman entitled ‘‘An Extension of 
the Theory of Resonances of Biological Cells: II. Cross-Section in 
a Plane Wave’? (Vol. 17, pp. 35-40). 

Equation (2) on page 36 should have a term 


substained from the right hand side, 

In equation (7), page 38, the numerical coefficient of T should be 
8 instead of 12, and on page 39, 24 instead of 36. These alter the 
numerical cross-sections but not their magnitude. Only the latter 
are quoted. 

Other minor errors on page 37 below equation (6): 
6 = II/2 should be 0, = 7/2. And on page 38 the next to the last 
line should read ‘‘lowest’? mode instead of ‘‘latest.” 
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Applications for the year 1957-1958 should be submitted before 
February 1, 1957, to the Fellowship Awards Board which is author- 
ized to make the seventh Sigma Delta Epsilon Fellowship Award of 
$1600. 

Women with the equivalent of a Ph.D. degree, carrying on re- 
search in the mathematical, physical, or biological sciences, who 
need financial assistance and give evidence of high ability and 
promise are eligible. During the term of her appointment the ap- 
pointee must devote the major part of her time to the approved re- 
search project, and not engage in other work for remuneration 
(unless such work shall have received the written approval of the 
Board before the award of the fellowship). 

Application blanks may be obtained from Dr. Dorothy Quiggle, 
Petroleum Refining Laboratory, The Pennsylvania State University, 
University Park, Pennsylvania, Announcement of the award will 
be made early in March 1957. 
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